1/2 


R0-IU.<4  3M 
UNCLASSIFIED 


DETECTION  OF  RTNOSPHERIC  CARBON  DIOXIDE  FRON  A 
SHUTTLE-BORNE  LIDAR<U>  AIR  FORCE  INST  OF  TECH 
HRIQHT-PATTERSON  AFB  OH  SCHOOL  OF  ENGINEERING  R  H 
DEC  02  RFIT/GE/PH/82D-2S  F/G  4/1  NL 


SELECTEj 
FEB  1  4  1960 

DETECTION  OP  ATMOSPHERIC 
CARBON  DIOXIDE  FROM  A 
SHUTTLE-BORNE  LIDAR 


THESIS 


Robert  B.  Wank 
AFIT/GEP/PH/82D-25  Cpt.  USA 


AFIT/GEP /PH/8  2  D- 2  5 


DETECTION  OF  ATMOSPHERIC 
CARBON  DIOXIDE  FROM  A 
SHUTTLE-BORNE  LIDAR 


THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 

Master  of  Science  - - 


Robert  H.  Wankr  B.S. 
Cpt.  USA 

Graduate  Physics 
December  1982 


Accesion  For 


NTIS  CRA&I 

—3— 

OTIC  TAB 

□ 

Unannounced 

□  1 

Justification 

- . J 

By . 

Dist  ibutioi'./ 

Availability  Codes 

I  Avail  and/or 
Special 


Approved  for  public  release;  distribution  unlimited 


PREFACE 


This  thesis  is  concerned  with  investigating  the 
detection  of  atmospheric  carbon  dioxide  at  various 
temperatures  and  altitudes  from  a  shuttle-borne  differential 
absorption  lidar.  Absorption  is  from  the  (ll'0-03'0)j 

CO2  band,  and  closely  parallels  the  earlier  frequency 
recommendations  of  Dr.  C.L.  Korb  of  the  NASA/Goddard  Space 
Flight  Center. 

I  would  like  to  thank  my  advisor,  Ma j .  James  J. 
Lange,  for  his  patience,  interest,  and  direction  in  the 
completion  of  this  thesis.  I  would  like  to  thank  Dr.  W.B. 
Roh,  Dr.  E.A.  Oorko,  and  Dr.  R.L.  Hengehold  for  their  input 
in  resolving  various  problems  that  I  encountered. 


L  I 


‘  •' ‘-'l 


Robert  He  Wank 


Contents 


Preface 


List  of  Figures 
List  of  Tables  . 


Abstract 


I.  Introduction 


Problem 


General  Approach 


Assumptions 


II.  Atmospheric  Carbon  Dioxide 

III.  Structure  of  Carbon  Dioxide 


IV.  Lidar  Detection  of  Carbon  Dioxide 


Lidar  Concept 


Lidar  Equation  Computations 
Mie  and  Rayleigh  Scattering 


Carbon  Dioxide  "Off "-Resonance 
Cross  Section  . 


Carbon  Monoxide  Cross  Section 


Water  Vapor  Cross  Section 


Ozone  Cross  Section 


Extinction  Coefficient  Calculations 


.  /.  /.  /.  1.S  .’ik*.  jiV 


Page 


Extinction  Due  to  Atmospheric  Absorption  .  . 

Carbon  Dioxide  "On "-Resonance 

Cross  Section . . . 


. ...t 


Extinction  Coefficient  of  "Off"- 

Resonance  CO2  . 

Extinction  Coefficient  of  Carbon  Monoxide 

Extinction  Coefficient  of  Ozone  .... 

Extinction  Coefficient  of  Water  Vapor 

Calculation  of  "Of f "-Resonance 

Return  Signal  Power  . 

Background  Noise  Calculations  . 

Signal-to-Noise  Calculations 

(Direct  Detection)  ...  . 

Heterodyne  Detection  . 

Detection  of  Hot  CO2  in  the  Atmosphere  . 

V.  Error  Analysis  of  DIAL  Calculations  .  .  . 

VI.  Conclusions  and  Recommendations  . 

Conclusions  .  . 

Recommendations  . 

Bibliography  . 

Appendix  A:  Rayleigh  Scattering  Coefficients 

Appendix  B:  Mie  Scattering  and  Extinction 

Coefficients  (4.86m)  . 

Appendix  C:  CO2  Absorption  Cross  Section 

^ ( "On" ) -Resonance  . 

Appendix  D:  CO2  Absorption  Cross  Section 

^ ( "Of f "-Resonance)  . 

Appendix  E:  CO  Absorption  Cross  Section 

( "Of f "-Resonance)  ....... 


Appendix  F; 


H.O  Absorption  Cross  Section 
( "Off "-Resonance)  .  .  .  . 


Appendix  G: 


0^  Absorption  Cross  Section 
( "Off "-Resonance)  .  .  . 


Appendix  H: 


COj  Extinction  Coefficient 
^ ( "Off "-Resonance)  .  .  . 


Appendix  I; 


CO  Extinction  Coefficient 
( "Off "-Resonance)  .  . 


Appendix  I; 


0^  Extinction  Coefficient 
( "Of f "-Resonance)  .  . 


Appendix  K: 


HyO  Extinction  Coefficient 
{ "Of f "-Resonance)  .  .  . 


Appendix  L: 


Hot  CO^  Absorption  Cross  Section 
{ "Or  f "-Resonance ) . 


Vita 


Coriolis  Forces  in  Linear  CO^  . 

Classical  Motion  of  Nuclei  in  Linear  CO2 

Simplified  Diagreun  of  the  Differential 

Absorption  and  Scattering  LIDAR  Technique 

Block  Diagram  of  a  Laser  Radar  With 

Heterodyne  System  . 


Abstract 


Remote  sensing  of  atmospheric  carbon  dioxide  from  a 
shuttle-borne  differential  absorption  lidar  (DIAL)  was 
investigated.  The  investigation  followed  the  basic  fre- 
quency  recommendations  of. Korb,  which  included  "on"- 
resonance  frequency  at  the  P(34)  line  of  the  (11 '0-03 '0)’ 

X 

absorption  band  of  C02.  Extinction  coefficients  for 
C02»  H2O,  CO,  Oj,  Rayleigh,  and  Mie  scattering  were 
determined  to  compute  estimated  return  signal  strengths  from 
various  altitudes.  Direct  detection  was  found  to  be 
unsuitable,  but  heterodyne  detection  provided  adequate 
signal-to-noise  ratios.  Cross  sections  for  hot  CO2  were 
also  calculated,  and  maximum  detectable  CO2  concentrations 
were  determined.  Finally,  an  error  analysis  of  predicted 
CO2  concentrations  was  made,  which  determined  that  due  to 
the  large  CO2  "on"-resonance  extinction  coefficient, 
accurate  measurements  of  CO2  below  five  kilometers  could 

nO^  b€  in£lCl0«  •  r  /  if  >  T' ■  * 
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Korbf  C.L.  et.  al,  "A  Lidar  Technique  for 
Measurement  of  Atmospheric  Carbon  Dioxide,"  Proceedings 
^he  Tenth  International  Laser  Radar  Conference. 
irrs-ltJel  sTlver  Springs,  Maryland:  Committee  on  Laser 
Atmospheric  Studies,  October  1980. 


viii 


DETECTION  OF  ATMOSPHERIC  CARBON 
DIOXIDE  FROM  A  SHUTTLE-BORNE  LIDAR 

I .  Introduction 

Atmospheric  carbon  dioxide  (CO2)  has  been  gaining 
significantly  in  abundance  at  various  monitored  locations 
since  the  turn  of  the  century  (Ref.  6:2).  Only  recently, 
however,  have  scientists  questioned  the  long-range  environ¬ 
mental  impact  of  this  increase.  Data  is  not  currently 
available  that  identifies  global  abundance  trends,  major 
sources  and  sinks  of  CO^,  and  transport  characteristics  at 
different  altitudes.  In  1979,  the  NASA  Atmospheric  Lidar 
Working  Group  submitted  a  report  (Ref.  7)  detailing  26  candi 
date  experiments  for  shuttle-borne  lidar  system,  one  of 
which  recommended  the  measurement  of  the  CO2  mixing  ratio  at 
different  altitudes. 

Little  actual  measurement  of  CO2  has  been  accomp¬ 
lished  using  a  long-range  lidar  system.  Remsberg  and 
Go'"^ley  (Ref.  15)  have  analyzed  the  detection  of  several 
atmospheric  constituents  in  both  the  troposphere  and  strato¬ 
sphere  from  a  shuttle-borne  system,  predicting  only  water 
vapor  profiles  as  attainable  in  the  troposphere.  Korb 
(Ref.  17)  has  proposed  an  aircraft  lidar  operating  at  4.86p 
using  a  frequency-doubled  CO2  laser  to  measure  CO2 
concentrations  in  the  lower  atmosphere. 


Problem 


The  problem  of  this  thesis  is  to  investigate 
different  lidar  systems  as  possible  candidates  of  a 
shuttle-borne  system  designed  to  measure  CO2  concentra¬ 
tions  at  various  altitudes  and  temperatures.  The  analysis 
is  directed  toward  applying  the  saune  general  design  para¬ 
meters  suggested  by  Korb  in  his  presentation. 


General  Approach 

The  investigation  centers  on  possible  lidar  systems 
capable  of  providing  data  from  a  shuttle-borne  system. 

Using  typical  shuttle  lidar  power  and  optics  parauneters,  and 
the  suggested  lidar  frequencies  of  Korb,  returned  signal 
powers  are  calculated  for  various  altitudes.  Signal-to- 
noise  calculations  are  then  made  for  each  altitude.  Hot 
CO2  concentrations  are  introduced  below  two  kilometers, 
and  a  maximum  detectable  CO2  concentration  is  determined. 
Finally,  an  error  analysis  of  the  predicted  CO2  concen¬ 
trations  is  completed. 


Assumptions 


The  atmosphere  is  considered  to  be  nonturbulent 
without  cloud  cover.  Visibility  is  considered  to  be  23 
kilometers  or  a  light  haze. 


Water  vapor  and  ozone  concentrations  are  taken  from 


the  mid-latitude  model  of  the  U.S»  Standard  Atmosphere, 

1976  (Ref.  39). 

Aerosol  concentrations  are  derived  from  the  model  by 
Shettle  and  Fenn  (Ref.  25). 


II.  Atmospheric  Carbon  Dioxide  Effects 


Carbon  dioxide  (CO2)  constitutes  approximately 
330ppm  of  the  atmosphere  by  volume  (Ref.  1:14).  Along  with 
CO,  1^2^'  ^^4'  ^2'  generally  considered  a 

uniformly  mixed  gas,  or  one  that  closely  maintains  that 
330ppm  concentration  over  vertical  profile  up  to  about 
eightly  kilometers.  Junge  and  Hagemann  are  commonly  cited 
for  their  studies  in  this  area  (Ref.  2:7).  However, 

Manquero  (Ref.  3:3)  reported  in  September  1975  variations  of 


a  CO2  vertical  profile  up  to  lOppm,  and  Hinkley  (Ref.  4:15) 
cites  CO2  concentrations  in  urban  areas  as  high  as  550ppm 
close  to  the  earth's  surface. 

Though  CO2  is  not  a  toxic  gas,  it  is  the  only  gas 
or  liquid  by-product  of  burning  fossil  fuels  that  has  been 
been  retained  in  large  quantities  in  the  atmosphere  (Ref.  5:9) 


In  fact,  experts  have  estimated  that  CO2  concentrations 


have  increased  from  fifteen  to  twenty  percent  since  the  turn 
of  the  century  (Ref.  6:2),  and  production  of  CO2  continues 
to  increase  at  a  4.3  percent  annual  rate  (Ref.  5:17).  The 
long  range  effects  of  increased  CO2  in  the  atmosphere  are 
purely  speculative.  Numerous  models  have  been  developed 
concerning  the  climatic  consequences  of  changes  in  the 
temperature  patterns  around  the  world.  The  most  common 
theory  suggests  that  an  increased  CO2  concentration  will 
result  in  a  general  warming  of  the  lower  atmosphere  and  a 
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cooling  of  the  stratosphere.  This  will  occur  due  to  absorp¬ 
tion  of  the  earth's  infrared  radiation  by  €©2#  and  subse¬ 
quent  emission  of  a  portion  back  into  the  atmosphere.  This 
radiation  trap  is  known  as  the  "greenhouse  effect"  (Ref.  6:38) 

One  climatic  model  predicts  doubling  the  CO2 
concentration  will  create  a  global  temperature  increase  of 
1.5®  to  4®C  by  the  first  half  of  the  next  century,  an 
increase  considerably  greater  than  any  resulting  from 
natural  causes  over  the  past  several  thousand  years 
(Ref.  6:45).  Another  model  predicts  a  temperature  increase 
from  1"C  by  2020  to  and  11®C  increase  by  2075  (Ref.  5:23). 

Kellogg  (Ref.  6:4)  has  listed  several  possible 
consequences  of  a  global  temperature  increase: 

1.  Middle  latitude  communities  will  require  less 
energy  demand  in  the  winter,  but  more  in  the 
summer . 

2.  A  1®C  change  could  add  ten  days  to  the  growing 
season,  and  increased  CO2  could  enhance  the 
photosynthesis  process.  However,  pest  losses, 
already  at  25  percent,  would  likely  increase. 

3.  Natural  ecosystems  would  shift,  as  some  species 
would  die  out,  and  others  migrate  to  new  areas. 

4.  Precipitation  patterns  would  shift,  decreasing 
rainfall  in  the  Midwest  of  the  U.S.,  and  the 
U.S.S.R.,  while  aiding  several  semi -arid 
developing  countries. 


5.  Climatic  changes  could  shift  ocean  currents,  sea 
surface  temperatures,  and  wind  patterns, 
altering  the  fishing  industry  production. 

6.  Labor  production  would  decrease  2  to  4  percent 
for  every  1®C  over  20®C,  and  diseases 
previously  confined  to  the  tropics  could  spread 
to  temperate  regions. 

7.  Melting  of  the  polar  ice  caps,  and  the 
subsequent  rise  of  the  oceanic  levels  could 
severely  affect  tourism,  and  result  in  major 
population  shifts.  For  instance,  a  five  meter 
rise  would  affect  eleven  million  people  living 
on  the  Atlantic  and  Pacific  coasts. 

Whether  or  not  any  of  these  predictions  will 
actually  occur  is  not  known.  Climatic  models  have  not  been 
developed  that  are  sophisticated  enough  to  consider  all 
possible  long  and  short  range  factors,  due  mainly  to  the 
lack  of  adequate  atmospheric  data.  Predicting  changes  to 
the  atmospheric  CO2  concentration  requires  investigation 
of  both  sources  and  possible  sinks.  Major  sources  of  CO2 
emission  include  fossil-burning  plants,  cement  manufacturing 
facilities,  and  possibly  major  CO2  laser  facilities.  The 
CO2  transport  process  is  not  yet  fully  understood.  Of  all 
the  estimated  CO2  produced  by  fossil-fuel  burning  during  a 
given  time  period,  about  50  percent  may  actually  be  retained 


in  the  atmosphere  (Ref.  7:11).  Others  estimate  that  the 


oceans  previously  absorbed  up  to  45  percent  of  all  CO2 
emissions,  and  forests  almost  all  of  the  remainder. 

However,  deforestation  has  resulted  in  the  oceans  being 
required  to  serve  as  a  sink  for  up  to  80  percent  of  the 
CO2  emissions,  a  percentage  they  are  unable  to  adequately 
absorb  (Ref.  6:37). 

Developing  an  accurate  climatic  model  is  necessary 
to  determine  whether  CO2  control  measures  should  be 
imposed.  Any  control  should  be  carefully  evaluated  since 
costs  will  be  extensive.  For  example,  it  has  been  estimated 
that  to  remove  50  percent  of  CO2  emissions  from  a 
conventional  coal-fired  plant  would  reduce  the  power 
generation  efficiency  from  34  percent  to  25  percent,  and 
double  the  cost.  Removing  90  percent  of  CO2  emissions 
would  decrease  efficiency  to  6  to  15  percent,  and  increase 
the  cost  by  a  factor  of  four  to  seven  percent  (Ref.  8:iii). 
However,  a  sophisticated  climatic  model  can  only  be 
developed  with  accurate  knowledge  of  CO2  sources,  sinks, 
and  transport  processes. 
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III.  structured  of  Carbon  Dioxide 

Carbon  dioxide  (CO2)  is  a  linear  triatomic 
molecule.  The  electronic  absorption  spectrum  is  in  the  far 
ultra-violet,  unavailable  for  atmospheric  study,  and  no  pure 
rotational  spectra  exists  due  to  a  lack  of  dipole  moment. 

For  a  triatomic  linear  molecule  normal  motions  include  three 
translational  modes,  two  rotational  modes  (degenerate),  and 
four  vibrational  modes.  The  available  energy  bands  are  a 
function  of  three  vibrational  quantum  numbers  describing  the 
four  vibrational  modes: 

1.  ~  symmetric  stretch  motion. 

2.  ^2  -  bending  motion  (degenerate). 

3.  ^^3  "  asymmetric  stretch  motion. 

However,  the  symmetric  stretch  vibration  (N)^)'  optically 
inactive  because  the  symmetry  of  its  dipole  moment  is 
constant  during  vibration.  The  main  '^2  ^3 

vibrational-rotational  bands  are  located  at  approximately 
667.40cm”^  (15p)  and  2349. 16cm”^  (4.3p),  respectively. 
Overtones,  hot  bands,  and  other  frequency  combinations 
result  in  additional  bands  centered  at  10.4,  9.4,  5.2,  4.3, 
2.7,  2.0,  1.6,  1.4  microns,  and  several  weak  bands  from  .78 
to  1.24  microns  (Ref.  4:36). 

However,  it  is  insufficient  to  describe  the  energy 
levels  based  solely  on  three  vibrational  quantum  numbers, 
and  the  rotational  quantum  number  (J).  Two  additional 


quantum  numbers,  1  and  r,  are  necessary  to  account  for 
energy  level  shifts  due  to  the  molecule's  rotational 
direction  with  respect  to  the  vibrational  bending  motion, 
Coriolis  forces,  and  Fermi  resonance.  The  quantum  number  1 
accounts  for  the  contribution  of  the  bending  mode  to  the 
angular  rotation,  and  has  values  of  1  *  V2,  ^2-2,  ...  0. 

If  a  molecule  does  not  rotate,  it  vibrates  with  exactly  the 
same  frequency  for  the  two  bending  modes.  However,  rotation 
about  an  axis  affects  the  moments  of  inertia  of  the  tvro 
modes  differently  so  the  bending  frequencies  are  no  longer 
equivalent.  In  addition  a  fictitious,  but  necessary  force 
must  be  introduced  in  order  to  study  a  rotating  coordinate 
system  while  disregarding  the  rotation.  This  is  the 
Coriolis  force  and  is  equivalent  to  F  »  2V  X  W,  where  W 
is  the  vector  angular  velocity  of  rotation  of  the  coordinate 
system,  and  V  is  the  vector  velocity  of  motion 
through  the  coordinate  system  (Ref.  9:30).  Vibrational 
frequencies  and  are  excited  by  Coriolis  forces 
when  the  bending  vibration  is  perpendicular  to  the  angular 
momentum  J.  If  the  bending  vibration  is  parallel  to 
J,  no  additional  vibrational  modes  are  excited  since 
the  Coriolis  force  is  zero  (See  Figures  1  and  2).  Because 
of  these  two  effects  the  two  degenerate  bending  energy 
levels  are  slightly  split,  which  is  called  1-type  doubling. 
The  two  sets  of  levels  are  designated  c  and  d  having 


different  rotational  constants.  When  1*1,  the  splitting 
is  the  most  important.  McClatchey  (Ref.  10:11)  states  that: 


When  the  linear  molecule  has  a  center  of  symmetry 
(as  CO-)  with  a  '*0  at  both  ends,  the  paired  atoms 
zero  nuclear  spin  cause  zero  statistical  levels  of  a 
given  parity. 

This  means  that  for  the  ground  vibrational  level  and  all 
other  levels  with  1  *  0  and  v^  =  even  (  II g"*”  symmetry), 
only  even-numbered  J  values  are  found.  When  1*0  and 
Vj  =  odd  (  I  symmetry),  only  odd  J  values 
exist.  If  1  is  greater  than  zero,  the  c  and  d  levels 
have  different  symmetry.  For  example  if  1  *  1  (  n^)  the 
c  sublevel  has  odd  levels  of  J,  and  d  has  even. 

The  rotational  constant  (Ref.  9:25)  for  a  linear 
molecule  can  be  evaluated  from 


B 


B. 


J(J+1)D 


(1) 


where 

B^  *  equilibrium  value  of  the  rotational  constant 

oc  *  change  in  B  due  to  excitation  of  i^^ 
vibration 

D  *  change  in  B  due  the  centrifugal  stretching. 
€ 


Normally,  oc^  can  be  determined  by  the  separation  of 
lines  of  adjacent  vibrational  states.  However,  perturba¬ 
tions  between  vibrational  states  invalidate  this  simple 
relationship.  Fermi  initially  explained  that  these 


2093.350,  respectively 


IV.  Lidar  Detection  of  Carbon  Dioxide 
Lidar  Concept 

In  September  1977  a  working  group  was  convened  by 
NASA  to,  among  other  goals,  identify  areas  of  theoretical 
and  laboratory  research  necessary  to  support  the  planned 
development  of  a  space  shuttle-borne  atmospheric  lidar 
(Ref.  7:iii).  One  of  the  areas  suggested  for  investigation 
was  the  CO2  mixing  ratio  in  the  troposphere. 

A  lidar  (light  detection  and  ranging)  operates  very 
similar  to  a  radar,  except  wavelengths  vary  from  the  near  UV 
to  the  middle  IR.  A  lidar  system  is  commonly  composed  of  a 
laser  operating  on  a  specific  wavelength,  whose  beam  propa¬ 
gates  through  the  atmosphere,  and  is  subsequently  recaptured 
by  a  detection  system.  A  lidar  system  can  have  the  detector 
located  with  the  laser  (single-ended  system)  or  at  another 
location  (double-ended).  A  spaceborne  system  would  be 
single-ended.  By  analyzing  the  returned  laser  signal, 
numerous  atmospheric  properties  can  be  determined.  These 
include  aerosol  and  thin  cloud  detection,  temperature 
profiles,  water  vapor  concentrations,  trace  gas  detection, 
and  wind  velocity  profiles.  In  addition,  spaceborne  lidar 
has  the  ability  to  obtain  global-scale  data  of  high  spatial 
resolution  in  a  short  time  period.  Pulsed  lidar,  has  an 
inherent  range  resolution  capability,  enabling  three- 
dimensional  depiction  of  an  atmospheric  target,  a  feature, 
very  difficult  to  accurately  obtain  from  a  passive  system. 


Lidar  systems  are  not  new.  Numerous  types  exist 
that  have  been  successfully  built  and  tested.  Types  can  be 
divided  into  two  categories,  direct  and  indirect  scattering. 
Direct  scattering  includes  a  basic  lidar  to  measure  aerosol 
scattering  in  the  atmosphere,  Raman  lidar,  and  resonance 
( fluoresence)  lidar.  Indirect  scattering  includes  differ¬ 
ential  absorption  lidar  (DIAL)  or  differential  absorption 
scattering  lidar  (DAS),  and  additional  types  of  Raman  lidar. 

Raman  lidar  involves  the  detection  of  molecular  scat 

tered  Raman-shifted  wavelengths.  The  scattered  radiation  is 

unique  to  the  particular  molecule  being  studied,  and  the 

intensity  is  directly  proportional  to  the  concentration. 

Raman  lidar  was  first  demonstrated  by  Kobayasi  and  Inaba 

(Ref.  12:139)  to  detect  SOj  and  CO2  from  an  oil  smoke 

plume  using  a  Q-switched  ruby  laser.  However,  a  Raman 

signal  is  directly  proportional  to  the  Raman  backscatter 
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cross  section,  which  typically  is  about  10  cm  sr 
Therefore,  it  is  not  useful  at  ranges  typical  from  a  space 
platform  to  the  troposphere. 

Resonance  lidar  utilizes  a  pulsed  laser  tuned  on  an 
absorption  line  of  the  molecule  being  investigated  to  excite 
the  molecule  into  a  higher  energy  state.  The  molecule  sub¬ 
sequently  fluoresces  with  a  known  quantum  yield  at  a  parti¬ 
cular  wavelength,  which  is  detected  and  analyzed.  The  cross 
section  for  a  resonance  interaction  can  be  several  magni¬ 
tudes  greater  than  a  Raman  cross  section,  making  it  suitable 


for  spaceborne  lidar  application.  Mcllrath  has  suggested 
the  following  species  possibly  suitable  for  investigation  by 
fluorescence  lidar:  Li,  Na,  K,  Ca+,  Fe,  Cl,  0,  OH,  NO,  BaO, 


and  NO2V  (Ref.  13:495).  In  addition,  Kweder  has  proposed 
HF  and  I2  as  promising  candidates  for  study  (Ref.  14:81). 

The  only  real  disadvantage  of  this  techniqu'^  is  the  problem 
of  quenching.  Since  an  excited  state  has  a  definite 
lifetime,  non-radiative  collisions  between  that  molecule  and 
additional  atmosphere  constituents  decrease  the  number  of 
actual  radiative  decays,  which  affect  the  sensitivity  and 
accuracy  of  the  return  signal  and  subsequent  concentration 
measurement.  In  fact,  Remsburg  and  Gordley  in  their  space 
shuttle  lidar  analysis  state  that  the  fluorescence  technique 
would  be  valuable  only  at  altitudes  in  the  upper  strato¬ 
sphere  (Ref.  15:624),  a  conclusion  also  reached  by  Mcllrath 
in  his  investigation  of  OH  (Ref.  13:497). 

The  only  technique  considered  feasible  for  tropo¬ 
spheric  study  of  molecules  is  the  differential  absorption 
lidar  (DIAL).  This  lidar  uses  a  tunable  laser  or  two  separ¬ 
ate  lasers  that  provide  electromagnetic  radiation  at  two 
distinct  wavelengths,  one  on  a  resonance  absorption  line, 
and  the  other  adjacent,  but  tuned  off  that  absorption  line. 
Using  a  remote  reflection  such  as  the  earth,  the  two  signals 
are  returned  to  a  detector  and  analyzed.  By  comparing  the 
two  signal  intensities,  an  average  concentration  over  the 


path  can  be  determined.  If  molecules  or  aerosols  serve  as 
scattering  surfaces  for  the  two  laser  beams,  a  vertical 
profile  of  the  pollutant  being  studied  can  be  made.  This 
three-dimensional  technique  is  called  range-resolved  differ¬ 
ential  absorption  lidar.  This  technique  would  appear  to  be 
the  only  suitable  method  available  to  measure  the  vertical 
distribution  of  CO2/  and,  therefore,  will  be  used  during 
this  study. 

Several  studies  of  possible  DIAL  experiments  from  a 
space  platform  have  been  made.  In  fact,  Remsburg  and 
Gordley  predict  the  only  important  gaseous  constituent  in 
the  troposphere  able  to  be  accurately  measured  is  the  water 
vapor  (Ref.  15:629).  However,  this  study  will  investigate 
CO2  concentrations  using  the  4.82pim  band,  which  has  been 
successfully  used  on  two  separate  investigations.  First, 
the  Air  Force  Geophysics  Laboratory  NKC-135A  Infrared  Flying 
Laboratory  conducted  detailed  spectral  solar  transmissivity 
studies  in  1970  over  Johnston  Island  at  various  altitudes  in 
the  4.0-5.3(x  region  demonstrating  the  4.86|i  band  was 
clearly  distinguishable  over  absorption  by  other  atmospheric 
constituents  (Ref.  16:2454).  In  1980  Korb  proposed  using  a 
frequency-doubled  CO2  laser  to  investigate  CO2  concen¬ 
trations  from  an  aircraft  platform  (Ref.  17:105).  This  report 


parallels  closely  the  basic  physical  parameters  set  forth  in 
the  Korb  report.  These  include  the  following  facts; 


1.  The  "on”  absorption  line  in  the  (11 '0,  03 *0)^- 
(00 °0)  band  is  the  P(34)  line  located  at 

2050 . 568cm~^.  To  operate  on  this  line,  a 
frequency-doubled  pulsed  CO^  TEA  laser  tuned 
to  the  P(42)  line  (2050 . 596cm~^ )  of  the 
(00“l)-(02®0)  band  can  be  used. 

2.  The  "off"  absorption  line  is  provided  by  a 

frequency-doubled  pulsed  isotope 

TEA  laser  operating  on  the  (R14)  line  at 
2057 . 024cm”^.  This  particular  line  was 
selected,  because  absorption  by  interfering 
species  other  than  CO2  was  determined  to  be 
the  same  for  both  "on”  absorption  lines  and 
"off"  absorption  lines. 

3.  A  GdGeAs2  crystal  is  used  as  the  frequency 
doubling  material  with  an  efficiency  of  ten 
percent. 

Lidar  Equation  Computations 

The  basic  lidar  equation  has  been  derived,  and 
tested  in  numerous  articles  and  experiments.  Probably  the 
most  complete  derivation  has  been  done  by  Measures 
(Ref.  18:1092)  whose  derivation  has  the  flexibility  to  be 


applicable  to  Raman,  fluoresence,  and  DIAL  lidars.  The 
Ildar  equation  utilizes  the  Bouguer's  (Beer's)  Law  (Ref.  4:49) 


I(X)  =  lQ(X)exp[-Na  (X)LJ 


(3) 


where: 

2 

I  (X)  =  Initial  radiance  (watts/m  -  sr)  of  a 
o 

monochromatic  signal  of  assumed  plane  wave. 

2 

I(X)  =  Final  radiance  (watts/m  -  sr)  of  the 
signal  measured  at  a  distance,  L. 

3 

N  =  Concentration  (1/m  )  of  absorbing  species 
assumed  uniform  across  L. 
o(X)  =  Absorption  cross-section  of  the  species 
wavelength,  X. 

This  law  assumes  that  a(N))  is  independent  of  both  the 
incident  radiation  intensity  and  the  concentration  of 
absorbing  molecules,  N. 

This  report  will  use  a  combined  terminology  of 
Thompson  (Ref.  19:4),  Binkley  (Ref.  4:76),  and  the  European 
Space  Agency  (Ref.  20:57).  The  lidar  equation  for  DIAL 
measurements  is  as  follows: 


C  T 


(r) 


o  2 


8i . (r)A^Ke 


expp2 


Ki  (r)  +  N, 


rn 


(r)  a 


.  (r))d^ 


(4) 


Returned  power  signal  (watts)  from  range 
at  time  t. 

Subscript  denoting  the  wavelength  of  the 
laser  where  i  =  1  is  "on"  absorption  line 
and  i  =  2  is  "off". 

Subscript  denoting  vertical  range  cell 
under  investigation. 

Transmitted  power  at  time  t^  (watts). 
Velocity  of  light  in  vacuum. 

Integration  time  of  the  returned  power 
signal  (sec). 

Volumetric  backscattering  coefficient 

of  the  atmosphere  Ocm'^sr"^ ) . 

2 

Effective  receiver  area  (km  ). 

Range  from  the  space  shuttle  to  the 
range  cell  under  investigation. 

Extinction  coefficient  of  atmosphere 
minus  extinction  due  to  CO2  (km~^). 

Number  density  of  CO2  molecules 
(km"^) . 

Absorption  cross  section  of  CO2  (km  ). 
Transmittance  of  the  optics  of  the  return 
signal. 

Frequency-doubling  crystal  efficiency  and 
Transmittance  of  the  optics  for  the  out¬ 
going  signal. 


The  integration  time,  ,  is  selected  based  on  the  desired 

vertical  size  of  a  range  cell.  The  size  of  a  cell  is  given 

by  (See  Figure  3).  A  cell  size  of  Ikm  has  been 

selected  for  this  study.  The  initial  signal  power  is  deter- 
0 

mined  from  where  is  the  energy  in  joules  and 

the  t*  the  pulse  duration.  For  use  on  a  space  shuttle,  the 
NASA  Atmospheric  Lidar  Working  Group  suggested  a 
pulsed  COj  laser  with  energy  from  1-10  joules  at  l-50Hz 
with  a  150ns  pulse  (Ref.  7:172).  Brockman,  in  his  analysis 
of  0^,  H2O,  and  NH^  measurements  by  a  pulsed  CO2 
spaceborne  lidar  uses  an  energy  of  5  joules  with  1000  shots 
(Ref.  21:564).  Megie  and  Menzies  in  their  UV  and  IR  DIAL 
comparison  study  assume  a  pulsed  CO2  laser  with  an  energy 
of  2  joules,  lOOn  sec  pulse  duration  with  a  lOHz  repetition 
rate  (Ref.  22:1181).  This  study  will  use  Megie  and  Menzies' 
pareuneters.  Generally,  the  telescope  receiver  diameter  is 
assumed  to  be  lOOcm  for  a  space  shuttle  lidar  with  an 
optical  system  efficiency  of  30  percent.  The  frequency¬ 
doubling  crystal  efficiency  is  assumed  to  be  10  percent. 

The  term  exp  “2  /  *  {  a  j^)dr  is  the 

fractional  transmittance  of  energy  along  a  two-way  path 
between  the  lidar  and  the  range  cell.  The  total  extinction 
coefficient  minus  CO2  absorption  (C^)  can  be  divided 
into  three  components: 


RECEIVER  GATE  DURATION 


Ci3(X) 


Extinction  coefficient  due  to 


continuum  absorption. 

=  Extinction  coefficient  due  to 
Mie  and  Rayleigh  scattering, 
s  Extinction  coefficient  due  to 
selective  absorption,  which  is 
highly  wavelength  dependent. 


Equation  (4)  is  the  generally  accepted  form  for  the 
DIAL  equation.  However,  this  form  assumes  a  laser  of 
monochromatic  frequency,  when  in  fact,  the  laser  has  a 
definite  line  shape.  For  a  CO2  TEA  laser  operating  at  one 
atmosphere,  the  full  width  at  half  maximum  frequency  is 
approximately  4.2GHz  at  300®K,  or  .14cm”^  (Ref.  52:17) 
with  a  Lorentzian  line  shape.  Therefore,  a  more  correct 
statement  for  the  returned  signal  power  must  include  a 
convolution  of  the  absorption  and  laser  line  shapes,  given 


as  follows: 


P.j(r,X)  -  /  P)^ij(r,X)dX 


-f-A  xe 

/  P^)^(X)Bi^(r,X)  exp 

R  j  O 


[-2  +  N^,Q^(r)  o.(r,X))  dr]dX 

Even  considering  that  Bj^j(r,X)  is  not  highly 
wavelength  sensitive  in  the  Infrared,  the  remaining  integral 
is  very  difficult  to  evaluate  analytically  or  numerically. 
Therefore,  this  study  will  assume  that  the  additional  error 


introduced  by  neglecting  the  shape  of  the  laser  pulse  will 
not  significantly  alter  the  final  results.  This  has  been 
assumed  extensively  in  the  past  giving  mixed  results  when 
theoretical  extinction  coefficients  are  compared  to  actual 
experimental  results.  For  example,  several  studies  have 
been  done  on  DF  and  CO  laser  attenuation  at  various  atmos¬ 
pheric  pressures,  temperatures,  and  compositions.  These 
studies  around  5^  should  give  an  indication  on  whether  the 
attenuation  of  a  laser  in  that  region  for  this  study  can  be 
approximated  by  Equation  ( 4 ) . 

For  the  DF  laser,  theoretical  calculations  of  extinc' 
tion  coefficients  due  to  N^O,  CO,  0^,  N2,  and  aerosols 
were  initially  done  by  McClatchey  and  Selby  (Ref.  53:2)  in 
1972  using  a  simple  Lorentzian  line  shape  for  the  inter¬ 
fering  species  to  determine  the  attenuation.  Several 
studies  have  been  done  to  experimentally  confirm  or  refute 
these  calculations.  Mills,  in  his  very  complete  theoretical 
and  experimental  study  in  1976  (Ref.  54:108),  compares  his 
experimental  results  of  DF  laser  attenuation  by  and 
those  of  two  other  experiments  to  the  calculated  extinction 
coefficients  of  N2O  assuming  absorption  cross  sections 
with  a  simple  Lorentzian  line  shape.  His  results  compare 
very  favorably  (1  or  2  percent  discrepancy)  with  theoretical 
calculations.  While  the  other  two  studies  do  not  agree  as 
well,  extrapolation  of  the  data  is  considered  the  primary 
reason,  and  not  the  theoretical  technique.  Hills  also  cites 


studies  by  Myers  of  attenuation  of  the  DF  laser  by  C02» 
which  had  discrepancies  up  to  four  orders  of  magnitude  from 
the  theoretical,  and  assumed  to  be  caused  by  a  weak  CO^ 
band  not  included  in  computations  (Ref.  54:121).  Mills  also 
compares  his  experimental  attenuation  of  DF  laser  radiation 
due  to  H2O  to  the  combined  attenuations  of  calculated 
H2O  plus  continuum  H2O  absorption  calculations  of  Burch 
(Ref.  54:143).  The  simple  Lorentsian  line  shape  for  absorp¬ 
tion  again  gave  generally  good  results  with  a  maximum  of  40 
percent  error.  Finally,  Deaton  investigated  DF  laser  attenu¬ 
ation  by  methane  (Ref.  55:302).  Results  indicated  that 
theoretical  attenuation  predictions  were  up  to  two  orders  of 
magnitude  less  than  actual  experimental  results.  This  was 
explained  by  the  theoretical  calculations  considering  only 
eight  strong  methane  absorption  lines  between  2710.96  and 
2712.11cm“^,  when  actually  71  lines  lying  between  2700  and 
2720cm~^  all  contributed  to  the  total  absorption. 

For  attenuation  of  CO  laser  radiation,  McClatchey 
again  did  the  first  major  theoretical  calculations,  which  he 
based  on  a  monochromatic  CO  laser  using  a  simple  Lorentzian 
absorption  line  shape  (Ref.  56:7)  for  attenuation  by  H2O, 

^^2'  ^3'  ^^4'  aerosols.  This  was  followed  by  a  study 
by  Ford  (Ref.  57:37),  who  analyzed  the  parameters  of  CO  laser 
attenuation  particularly  effects  of  water  vapor  and  attenua¬ 
tion  at  the  extreme  wings  using  a  simple  Lorentzian  absorp¬ 
tion  line  shape.  This  study  was  continued,  comparing  Ford's 


experimental  results  of  CO  laser  attenuation  by  H2O  at  four 
different  pressures  to  calculated  attenuations  (Ref.  58:12). 
Generally,  results  were  quite  close,  with  any  discrepancies 
caused  by  failure  to  consider  absorbing  lines  or  effects  of 
extreme  wings.  Long  followed  up  this  study  with  additional 
results  in  1973  (Ref.  59:34).  His  data  showed  that  the 
attenuation  of  CO  laser  radiation  predicted  by  the  simple 
Lorentzian  line  shape  of  H2O  molecules  gave  consistently 
lower  values  than  experimental  values  of  extinction 
coefficients.  A  second  line  shape  was  derived  to  account 
for  the  additional  absorption  which  was  due  to  absorption  in 
the  extreme  wings. 

These  studies  of  CO  and  OF  laser  radiation  clearly 
indicate  that  even  though  the  laser  line  shape  has  not  been 
considered  in  the  past,  accurate  results  of  predicted  laser 
attenuation  can  be  made.  Any  discrepencies  between  experi¬ 
mental  results  and  theoretical  calculations  appear  to  be  due 
to  failure  to  consider  weakly  absorbing  lines,  absorption  in 
the  extreme  wings,  and  the  inadequacy  of  the  simple 
Lorentzian  line  shape  to  describe  extreme  wing  effects 
correctly.  However,  adding  a  factor  to  approximate  the  line 
shape  of  the  laser  can  only  confuse  the  situation,  and  may, 
in  fact,  lead  to  more  inaccurate  results.  Therefore,  this 
study  will  continue  to  use  the  convention  for  the  DIAL 
equation  given  in  Equation  (4). 


In  order  for  a  return  signal  to  be 

detectable  two  conditions  must  be  met  (Ref.  20:59): 


1.  P,  .  >  P  where  P  is  the  equivalent 
13  —  n  n 


noise  power  of  the  system. 


^2j  "  ^Ij  -  ^n 


If  it  assumed  that  ^  ^  ^  2  previously  stated  from 

the  Korb  study,  then  Equation  (4)  yields: 


where  Ac=  ^  1  ~  2 

Therefore,  the  first  condition  above  yields: 


P2j  exp 


"coj 


(5) 


(6) 


However,  the  power  signal-to-noise  ratio  is  defined  as  on  the 
"off"  absorption  wavelength: 


n 


S 

"n" 


(7) 


Combining  Equations  (6)  and  (7)  yields: 

Ac  d^ 


exp 


(8) 


The  second  condition  above  combined  with  Equation  (5)  gives: 


This  can  be  rewritten  into  the  more  useful  form: 


“p[-2  ScOj  i  ^  -  WSiTJ 

Combining  Equations  (11)  and  (8)  yields: 

^  "  {S/n)2  -  **002^^*^^  -  (S/N)2 


(11) 


(12) 


In  order  to  satisfy  both  conditions  of  detectability  it  can 
be  seen  that^— must  equal  two  as  a  minimum.  Therefore, 
it  is  necessary  to  evaluate  P2j  and  to  determine 
whether  the  parameters  of  this  particular  lidar  are 
acceptable  to  detect  atmospheric  C02« 

Mie  and  Rayleigh  Scattering 

Rayleigh  scattering  is  classical,  elastic  scattering 
involving  incident  electromagnetic  radiation  and  atmospheric 
molecules.  No  appreciable  loss  of  energy  occurs,  and  the 
frequency  is  not  changed.  Since  all  molecules  contribute  to 
Rayleigh  scattering,  conclusions  about  concentrations  of 
individual  constituents  cannot  be  made.  The  extinction 
coefficient  due  to  Rayleigh  scattering  is  given 
by  (Ref.  4:52): 


(13) 


CiR(X)  =  [8  -  l)^/3NgX^]  [(6  +  3^)/(6  -  7^)] 

where: 

N  =  Number  density  of  gas  molecules. 

g 

n  =  Refractive  index  of  the  medium, 
i  =  Depolarization  factor  of  scattered 
radiation  I  .035. 

X  =  Wavelength  of  incident  radiation. 
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Note  the  X  dependence.  Experiments  with  several 
atmospheric  models  has  yielded  a  simple  empirical  equation 
for  the  cross  section  due  to  Rayleigh  scattering  up  to  an 
altitude  of  100km  (Ref.  4:98): 


CiR(X)  =  4.56(X(H)/.55)“^  X  10"^^  cm^ 


(14 


Since  =  N^  u^r,  the  Rayleigh  extinction  coefficient 
can  be  calculated  for  various  altitudes,  using  number  densi 
ties  from  a  U.S.  Standard  Atmosphere  (Ref.  23:6).  See 
Appendix  A  for  computations.  After  calculating  the  extinc¬ 
tion  coefficient,  the  volumetric  backscattering  coefficient 
(8.  ._)  due  to  Rayleigh  scattering  can  be  computed.  8. 

1 jK  XJR 

is  a  measure  of  the  eunount  of  radiation  scattered  into  a 


solid  angle  from  a  given  volume  of  gas.  The  magnitude  of 
8. is  given  by  the  Rayleigh  scattering  law  (Ref.  20:24): 


Mie  scattering  requires  a  much  more  complex  model. 

It  involves  interaction  between  electromagnetic  radiation, 
and  particles  of  size  approximately  equal  to  or  larger  than 

the  wavelength  of  the  radiation.  Particle  size  can  range 
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from  10  p  for  the  dust  or  pollen  to  10  p  for  raindrops. 

Size  of  particles  particularly  affects  the  angular  distribu¬ 
tion  of  scattering  (Ref.  1:70).  The  interaction,  dominated 
by  elastic  scattering,  also  includes  an  absorption 
component.  An  additional  significant  factor,  besides  wave¬ 
length  and  particle  size  distribution,  is  refractive  index 
of  the  particle.  Due  to  the  complex  nature  of  Mie 
scattering,  the  general  approach  to  analyzing  its  effects  is 
the  use  of  an  empirical  atmospheric  model.  Several  exist 
including  Oeirmendjian  haze  models  (Ref.  26:187),  Elterman's 
Atmospheric  model  (Ref.  23:1),  a  simplified  LOWTRAN  MODEL 
(Ref.  27:11),  Shettle  and  Fenn's  model  of  atmospheric 
aerosols  (Ref.  25:2-1),  and  a  simple  empirical  equation 
relating  the  Mie  scattering  coefficient  to  visual  range: 

1' 

CiM  =  (3.91/V)  (.55/X)*^®^^’  km"^  (16) 

where  V  is  the  visual  range  in  kilometers  and  \  the 
wavelength  in  microns.  This  equation  is  often  used,  but 
would  be  difficult  to  analyze  for  various  altitudes,  and  the 
validity  of  its  extension  to  the  infrared  from  the  visual 
spectrum  has  been  questioned  (Ref.  24:1500).  For  this  study 
Shettle  and  Fenn's  aerosol  model  will  be  used  due  to  its 


extensive  experimental  basis  and  close  agreement  with 
Elterman's  model.  Penn's  model  does  not  include  turbulent 


effects  on  scattering  or  absorption.  The  study  divides  the 
atmosphere  into  four  regions: 

1.  A  boundary  layer  up  to  two  kilometers  describing 
aerosols  in  a  rural,  urban,  and  maritime 
environment. 

2.  Troposphere  up  to  about  ten  kilometers 
representing  spring  and  summer  or  fall  and 
winter  conditions. 

3.  Stratosphere  up  to  thirty  kilometers  for 
background,  moderate,  high,  and  extreme  volcanic 
conditions  for  two  seasonal  models. 

4.  Upper  atmosphere  including  a  model  with  back¬ 
ground  conditions,  and  one  with  a  high  aerosol 
concentration  for  those  altitudes. 

The  Shettle  and  Fenn  model  plots  the  attenuation  coefficient 
versus  wavelength.  By  combining  the  data  from  both  graphs, 
extinction  coefficients  and  scattering  coefficients  due  to 
Mie  scattering  can  be  estimated  at  4.82pm.  This  study  made 
the  following  assumptions: 

1.  Urban  model  in  the  boundary  layer  for  a  moder¬ 
ately  clear  atmosphere  (23km  visibility). 

2.  Only  normal  background  is  considered  in  the 
stratosphere,  and  volcanic  effects  are  neglected. 
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3.  Any  attenuation  over  30km  is  assumed  to  be 


negligible. 

Results  due  to  extinction  and  scattering  due  to  aerosols  can 
be  found  in  Appendix  B. 

Extinction  Due  to  Atmospheric  Absorption 

The  extinction  coefficient  due  to  absorption  can 
also  be  defined  as: 

Ci(X)  *  Ng  a.  (17) 

where  is  the  absorption  coefficient  for  a  molecule  at 
either  the  "on"  or  "off"  absorption  wavelength  and  is 
the  number  density  of  the  particular  molecule.  The  cross 
section  must  be  determined  for  three  different  cases. 

First,  in  the  lower  atmosphere,  collision  broadening 
dominates  and  the  absorption  line  generally  takes  the  form 
of  a  simple  Lorentzian  (Ref.  4:32): 

(SA  )ac 

a(\J)  =  - ^ - 5-  (18) 

(V  -  +  V 

where 

2 

a(\j)  s  Absorption  line  cross  section  (cm  ). 

S  s  Line  intensity  of  an  individual  spectral  line 
(cm^  cm“^  mol”^),  and  s  ■  /“  o(\))d\). 
ac_  >  Half  width  at  half  maximum  of  the 

li 

absorption  line  (cm”^). 

\)  s  Wavenumber  (cm~^)  of  incident  radiation. 

*  Wavenumber  (cm“^)  of  the  line  center. 


Quantum  mechanics  provides  this  expression  for  the  intensity 
of  separate  transitions  of  a  molecule  from  state  k  into 
state  1  (Ref.  4:34): 


®kl  “  (8iT\j^j^/3hc)  (N^/g^N) 


exp(-hc\) 


where 


=  Transition  wavenumber. 

Nj^/N  =  Fractional  number  density  in  state  k 
with  statistical  weight  gj^. 

=  Square  of  the  matrix  element  of  the 
dipole  moment. 

h  =  Planck's  constant, 

c  =  Speed  of  light, 

k  =  Boltzmann's  constant. 

T  s  Temperature  of  the  medium. 

The  simple  Lorentzian  is  not  the  only  line  shape  used  to 
describe  collisional  broadening.  Serious  problems  in  the 
extreme  wings  where  |\)  -  \)^|  is  very  large  has 
resulted  in  development  of  modifications  to  the  simple 
Lorentzian  such  as  the  Gross,  Van  Vleck-Weisskoph,  and  Van 
Vleck-Huber  line  shapes  (Ref.  28:9).  However,  none  of  these 
are  totally  satisfactory  for  extreme  wing  measurements. 
Therefore,  this  study  will  assume  a  collisional-broadening 
line  shape  of  a  simple  Lorentzian. 
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Collisional  broadening  dominates  below  about  20 
kilometers.  Above  30  kilometers,  however,  Doppler 
broadening,  which  is  independent  of  pressure,  dominates. 

The  absorption  cross  section  for  a  Doppler-broadened  line 
shape  is  Gaussian  given  by  (Ref.  4:33): 

o(\))  =  (S/aCjj)(ln  2/Tr)^  expj^-(\J  -  ln2/(Xj^^J  (20) 


where 


oCp  =  (\)Q/c)(2kT  In  2/m) 


where 


m  s  Mass  of  an  individual  molecule  of  a  particular 


It  will  be  shown  that  absorption  above  30  kilometers  is 
negligible  compared  to  effects  of  the  lower  atmosphere. 

For  regions  between  about  20  to  30  kilometers, 
absorption  line  shapes  are  a  convolution  of  the  (Gaussian  and 
Loren tzian  known  as  the  Voigt  line  shape.  The  Voigt  profile 
is  given  by  (Ref.  4:246): 


<7(^0 )  -  (S/A'0)Re  {w(x  »  iy)} 


where 


Re{w(x  +  iy}}  »  Real  part  of  the  error  function  of 
complex  argument,  w(x). 

A\)  *  (Trln  2)^ac. 


(N)  -  N)o)AN} 


Ol/aN) 


However,  to  simplify  calculations  in  this  study  a  simple 
Lorentzian  was  assumed,  which  is  a  valid  assumption  since 
almost  all  of  the  significant  absorption  occurs  below  10 
kilometers.  The  effect  of  using  a  Lorentzian  above  20 
kilometers,  if  anything,  results  in  a  more  conservative 
return  signal  than  if  the  Voigt  profile  was  included. 

The  total  extinction  coefficient  due  to  absorption 
is  given  by: 

m 

5i  -  H  Sik 


where 


m  a  Number  of  absorbing  species  in  the  atmosphere. 


The  Infrared  Handbook  (Ref.  29; 5-92)  provides  a  diagreun  of 
major  atmospheric  constituents  affecting  transmission  in  the 
4.86p  region.  These  include  CO^,  water  (H^O),  ozone 
(0^),  and  carbon  monoxide  (CO).  Effects  of  the  isotope 
(1-0)  band  are  neglected  due  to  its  small 
atmospheric  concentration.  For  this  same  reason,  effects  of 
the  (00°!  -  00^0)  band  of  carbonyl  sulfide  (OCS)  with  a 
band  center  at  2062. 2cm*'^  (Ref.  30:1325)  will  be  ignored. 

The  extinction  coefficient  for  the  major  absorbers  listed 
above  will  be  discussed  individually. 


Carbon  Dioxide  "On”  -  Resonance  Cross  Section 

The  basic  structure  of  carbon  dioxide  has  been 
previously  discussed.  A  major  effect  of  the  (ll'O  -  03*0)^ 
-  (00 ®0)  band  is  the  large  Coriolis  effect  resulting  from 
1-1.  These  Coriolis-type  resonances  induce  a  large  part 
of  the  transition  moment  in  the  R-  and  P-  branches,  thus 
"borrowing"  intensity  from  stronger  parallel  bands 
(Ref.  10:26).  Determination  of  actual  intensities  of 
individual  absorption  lines  has  not  always  been 
theoretically  predicted  or  experimentally  verified. 
Formally,  the  intensity  of  any  line  at  wavenumber  N)  can 
be  expressed  as  (Ref.  10:13) 


^  'Sv  ’SRot 


where 

O 

*  vibrational  intensity  of  a  nonrotating 
molecule  at  a  vibrational  origin, 

F  a  Factor  that  accounts  for  simultaneous  vibration 
and  rotation. 

For  a  first  order  theoretical  correction  the  F-factor  can  be 
expressed  as 


F  »  S 


nonrigid'^^rigid  * 


where 


m  «  (J"  +  1)  for  R-branch;  (-J")  for  P-branch;  (J") 
for  Q-branch. 


C=  Empirical  value  used  by  McClatchey  to  describe  a 


particular  vibrational  band. 

As  further  simplification  for  C02»  McClatchey  used  the 
equation  below  to  correct  for  the  Coriolis  effect,  and 
determine  individual  line  intensities: 

S„  =  S,  (1  +  m  C (25) 
m  V 

where  for  the  (ll'O  -  03 *0)^  -  (00“0)  band  (Ref.  31:518) 

“  -22  -1  2  -1 
=  35.9  X  10  ^mol  -^cm^cm 

C  =  -.041 

However,  this  expression  would  appear  to  give  erroneous  data 
for  the  P-branch  since  S_  would  continue  to  increase  as  J 

m 

increases.  Burch  (Ref.  32:4-7)  completed  excellent  high 
resolution  spectroscopy  of  this  band  in  1970,  which  confirms 
the  predicted  reduction  of  the  R-branch,  and  enhancement  of 
the  P-branch.  Snyder  (Ref.  32:3)  reported  that  an 
additional  higher  order  term  of  the  F-factor  must  be  used 
for  the  (ll'O  -  03 *0)^  -  (00*0)  band  so 

P  =  (1  +  Cj^m  +  C  (26) 

andCj^  »  -.0237,  ^2  “  -.000527 

This  analysis  gives  results  which  more  accurately  coincide 
with  the  spectra  of  Burch  than  the  simple  first  order 
approximation.  Experimental  verification  of  individual  line 
half-widths  has  not  been  completed.  Therefore,  this  study 


will  use  an  average  value  of  =*  .07cin  for  one 

atmosphere  of  pressure  at  296*‘K  (Ref.  10:6).  Using 

Equations  (18)  and  (26)  with  the  "on”  wavenumber  of 

2050.596cm'~^,  a  total  absorption  cross  section  was 
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computed:  ^2C0  *  2.029  x  10  cm  . 

Appendix  C  lists  the  individual  line  strengths  and 
absorption  cross  sections  that  significantly  contribute  to 
the  total  absorption.  Note  the  dominate  absorption  cross 
section  of  the  P(34)  line. 


Carbon  Dioxide  "Off”  -  Resonance  Cross  Section 


This  analysis  parallels  closely  to  that  of  the  "on 
wavelength  except  \)  =  2057 .024cm~^,  corresponding  to 
the  frequency-doubled  ^^C^^Oj  laser  (00®1  -  (10*0/  02*0 )j.j 
Again,  oc^  «  .07cm~^.  Using  Equations  (18)  and  (26)/ 
the  total  absorption  cross  section  can  be  computed. 
Appendix  D  lists  the  individual  line  intensities  and  cross 

sections  for  CO^  at  off  resonance  with  the  total  ^2CO  * 

-21  2  ^ 

X  10  cm  .  This  is  approximately  15  times  less  than  the 

on-resonance  cross  section. 

Carbon  Monoxide  Cross  Section 

Carbon  Monoxide  (CO)  is  a  well-studied  diatomic 
molecule.  In  general  for  a  diatomic  molecule,  molecular 
rotational  transitions  (  J  «  ^1)  have,  to  a  good 
approximation,  the  same  transition  probability.  However, 


absorption  lines  do  not  have  the  same  intensity.  Instead, 


intensities  are  directly  proportional  to  the  initial  number 
of  molecules  in  each  level  (Ref.  51:11-12).  The  (1-0)  band 
of  CO  with  a  band  center  at  2134.2716cm''^,  dominates  the 
spectrum  in  the  4.86p  region.  The  band  sum  intensity, 
is  reported  by  Rothman  (Ref.  34:793)  as  981.3  x  10~^®cm”^mol'’^cm 
In  order  to  determine  the  fraction  of  the  total  band 
intensity  that  each  individual  absorption  line  contributes, 
molecular  number  densities  for  each  rotational  level  of  the 


vibrational  ground  state  must  be  computed.  The  relationship 
between  .the  number  density  in  a  particular  vibrational  level 
to  the  total  number  density  assuming  a  Boltzman 
distribution,  is  given  by 


where 


\  =  he 


[  “e 


+  i)  -  “g  (V  +  i)^ 


Q^( vibrational  partition  function) 

N  -he  CO  (V  +  J)  -  (0  (V  +  i)^  Xo 

»Z  exp  - 2 -  - e - e 

v=o 


0)  «  Hypothetical  equilibrium  frequency  of  an 

anharmonic  oscillator  and  for  CO  has  a  value  of 
2169. 7 4cm"^  (Ref.  51:V-11) 


Xg  ■  Anharmonicity  constant  for  a  diatomic  molecule 
and  has  a  value  of  .0061  for  CO  (Ref.  51:111-20) 


To  compute  it  can  be  shown  that  essentially  all 


molecules  are  in  the  ground  vibrational  level  at  a 
temperature  of  296®K 
n. 


n 


-  exp/  3.086  X  lO"^ 


v=0 


v=2 


=  1.05  X  10 


kT 

-9 


v=0 


(30 


(31 


Like  carbon  dioxide,  carbon  monoxide  is  considered  by 

McClatchey  as  a  uniformly  mixed  gas  in  the  atmosphere,  with 

12  16 

a  number  density  of  .075ppm.  The  C  0  isotope 
composes  98.652  percent  of  the  atmosphere  CO  (Ref.  10:6). 
From  the  U.S.  Standard  Atmosphere  (Ref.  23:6)  the  total 
molecule  number  density  at  the  earth's  surface,  is 

2.547  X  IO^^CM"^.  Therefore, 

12pl6.  «3 

**Tot  “  *  lO'^’^cm 

All  of  the  molecules  must  occupy  a  particular 

vibrational  energy  level,  and  so 


N 


Tot 


n 


v»0 


*  “v-l  *  "v-2  * 


Combining  with  Equations  (30)  and  (31) 


-5  -9 

n^^^  ^  :  n^Q  (1  +  3.086  x  10  ^  +  1.05  x  10  ^) 


:  1.884  X  lO^^cm"^ 


n _  ”  n  . 


since  essentially  all  of  the  CO  molecules  are  in  the 
ground  vibrational  state,  it  is  now  necessary  to  determine 
their  distribution  in  the  rotational  levels.  Again, 
assuming  a  Boltzman  distribution  and  including  the 
degeneracy  due  to  angular  momentum,  the  molecular  number 
densicy  of  a  rotational  level  to  the  total  number  density  of 
that  vibrational  band  is  given  by 
n_  {(2J  +  l)exp(-E_AT) } 

J  ^  'J.  _  /^*i\ 


where 


(Rotational  partition  function) 

N 

+  l)exp(-Bhc  J(J  +  D/lcT)}  (34) 

J*0 

B  (Rotational  constant)  ^  1.9225cm~^  (Ref.  35:207). 


The  factor  of  i  is  needed  since  only  the  P>branch  of  the 
band  contributes  to  the  absorption.  must  be  computed 
for  each  value  of  J.  For  this  study  computed  values  for 
J  «  0  to  J  ■  30  were  used.  Values  below  J  *>  30  were 
insignificant.  Summing  over  the  31  values  gives  I  107.32 
Knowing  the  total  band  intensity  and  using  Equation  (33), 
the  individual  absorption  line  intensities  can  be  computed 
from  the  simple  proportion 
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Again,  since  individual  absorption  line  half-widths  are  not 

available,  an  average  value,  ot^^,  of  .OScm”^  (Ref.  10:6) 

is  used.  Using  Equation  (18)  where  'O  =  2057.024cm 

the  "off"  resonance  wavenumer,  individual  absorption  cross 

sections  can  be  computed.  Appendix  E  lists  these  results 

for  significantly  contributing  absorption  lines,  along  with 

the  line  centers  (Ref.  35:204),  and  line  intensities.  A 

summation  of  the  individual  line  cross  sections  provides  a 
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total  cross  section  of  ^ 2C0  *  3.585  x  10  cm  . 

Water  Vapor  Cross  Section 

Water  has  the  structure  of  an  asymmetric- top 
molecule,  which  is  a  rotor  where  no  two  principal  moments  of 
inertia  are  equal.  This  necessitates  the  use  of  three 
quantum  numbers  to  describe  a  particular  rotational  level 
(Jf  K  ,  K  )  (Ref.  10:12).  K  and  K  assume  all  values 

SC  A  c 

0,  1,  2,  .  .  .  J,  subject  to  K  +  K  =  J  or  J  +  1. 

A  O 

Therefore,  2J  +  1  levels  exist  only  for  a  given  j.  Also,  a 
nuclear-spin  statistical  weight  must  be  applied  to  all 
asymmetric  top  molecules  with  two  or  more  identical  nuclei. 
Therefore,  the  probability  of  occupation  of  a  particular 
rotational  level  must  be  multiplied  by  a  weighting  factor  of 
3  or  1  depending  on  the  spin  of  the  nuclei  in  relation  to 
the  rotation  of  the  molecule  (Ref.  9:103).  These  and  other 
complexities  hinder  analysis  of  absorption  spectra. 
Fortunately,  a  recent  work  by  Flaud  and  Camy-Peyret  (Ref.  36) 


by  the  quantum  numbers  J 


predicts  at  \)  =  2056.0cm  ^  only  .879  of  the  absorp¬ 


tion  predicted  by  MWW  at  a  temperature  of  308 “K.  However, 
the  computed  value  of  the  total  cross  section  is  probably 
accurate  to  about  a  factor  of  two,  which  will  not  have  any 
significant  effect  on  the  total  extinction  coefficient  due 
to  all  of  the  interfering  species. 

Ozone  Cross  Section 

Ozone  (0^)  is  also  an  asymmetric-top  molecule, 
very  similar  in  structure  to  water.  However,  the  nuclear- 
spin  statistical  weight  factor  for  isotopically  symmetrical 
0^  has  values  of  0  or  1  for  the  odd  or  even  parity  of  J  + 

K_  +  K  +  V,  (Ref.  10:12).  The  primary  bands  contribu- 
ting  to  absorption  in  the  4.86p  region  include  the  (10®1) 

-  (00®0)  band,  centered  at  2110.785cm~^  with  a  band  of 

r  -20  -1  -1  2 

Sjjj  =  113.4  X  10  cm  mol  cm  ,  and  the 

(00*’2)  -  (OO^O)  band,  centered  at  2057.892cm"‘^  with  a  band 
strength,  2  Sjjj  *  11.1  x  lO'^^cm-^mol’^cm^.  The  (20"0)-(00*0) 
band  centered  at  2201.157  also  contributes  slightly 
(Ref  34:792).  Again,  it  is  fortunate  that  Flaud  and  Czuny- 
Peyret  in  another  study  have  evaluated  this  particular 
molecule,  identifying  transitions  by  the  three  rotational 
quantum  numbers,  their  respective  wave  number,  and  the  inten¬ 
sity  of  the  transition  (Ref  37:192).  A  line  half-width, 
oc^  s  .llcm~^  is  assumed  (Ref  10:6).  Again,  using 
equation  (18)  with  \)  2057 .024cm~^,  individual  cross 

sections  can  be  calculated.  Appendix  G  lists  the  wave 
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number,  intensity,  and  cross  section  for  170  lines  that 

contribute  significantly  to  the  total  absorption  cross 

section.  Summing  individual  cross  sections  gives  a  total 
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cross  section  of  =  1.37  x  10  cm  . 


Extinction  Coefficient  Calculations 

Cross  sections  have  been  calculated  for  the 
interfering  species  of  C02»  CO,  H2O,  and  0^.  However, 
these  calculations  are  made  at  latm  and  296“K.  In  fact  the 
line  intensity  is  a  function  of  temperature.  Therefore,  if 
the  Lorentzian  line  shape  is  valid  up  to  30  kilometers  in 
altitude,  cross  sections  as  a  function  of  altitude  must  be 
calculated.  If  these  cross  sections  are  known,  then  using 
Equation  (17)  extinction  coefficients  can  be  calculated  by 
multiplying  the  cross  section  by  the  interfering  species' 
number  density  as  a  function  of  altitude. 

The  intensity  of  absorption  line  as  a  function  of 
temperature  is  given  by  McClatchey  (Ref.  10:3)  as: 


S(T) 


S(Tg)Q^(T3)Q^(Tg) 

Q^(T)Qj.(T) 


exp 


1.439E"(T  -  T„) 
s 

ttI 


(36) 


where 


E"  =  Energy  of  the  lower  state  of  the  transition 
(cm“^) . 

Q^,Q^  ®  Vibrational  and  rotational  partition  functions 
S(T  )  » 


Line  intensity  at  the  standard  temperature. 


Tuer  (Ref.  1:16)  has  reported  the  following  equation  rela¬ 
ting  the  line  half -width  to  temperature  and  pressure: 


aCj^(T,P) 


«T (P/P„) (T^/T)" 
Xj  s  s  s  s 


(37 


where 

ac_(T  P  )  =  Line  half-width  at  temperature, 
li  s  s 

T  =  296 “K,  and  pressure,  P  =  1  atm. 
s  s 

n  s  assuming  temperature  independent  collision 
diameters . 

This  value  for  n  is  uncertain  since  values  have  been 
reported  up  to  n  -  1.  However,  n  =  i  is  the  most  commonly 
adopted  figure,  and  will  be  used  in  this  study. 

From  Equation  (18)  it  can  be  seen  that  the 
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absorption  line  cross  section  varies  as  oc^  varies  in 
the  denominator.  However,  the  contribution  of  this  term  is 
insignificant  unless  the  incident  radiation  wavenumber 
happens  to  coincide  very  closely  with  an  absorption  line 
center,  as  is  the  case  when  operating  on  the  ”on”-resonance 
laser  frequency.  However,  when  operating  *off "-resonance, 
this  factor  can  be  ignored.  Therefore,  by  combining 
Equations  (18),  (35),  and  (36),  an  "off "-resonance  cross 
section  as  a  function  of  temperature  and  pressure  can  be 
calculated  from: 


fs 


Qv(Ts) 


2  Q„(T) 


QrCTs)"]  jl.439{T-Tg 
Q^(T)  1  TtI 


(38) 


Extinction  Coefficient  of  "Off"  -  Resonance  CO^ 

The  vibrational  band  of  the  lower  energy  level  of 
the  (ll'O  -  03'0)j  -  (00®0)  transition  is  ground  state. 

The  energy  for  a  particular  rotation-vibrational  level  for  a 
linear  triatomic  molecule  can  be  approximated  as  (Ref.  51: VI-6) 

I  B"J"(J"  +  1)  +]^  (v?  +  J)\)^  ( 

where  for  the  ground  vibrational  state 
B"  *  .3902cm"^  (Ref.  11:137) 

\)j^  =  1330. Ocm“^;  \)2  *  667.3cm"^j  *  2349. 3cm"^  (Ref.  33:VI-6) 

Since  it  would  be  very  difficult  to  calculate  an  individual 
line  absorption  cross  section  as  a  function  of  altitude  for 
each  line  contributing  to  absorption,  a  rotational  quantum 
number,  J",  will  be  selected  corresponding  to  the  line  that 
contributes  most  significantly  to  the  total  absorption.  In 
this  case  that  line  is  the  P(26)  line  (See  Appendix  D) ,  so 
J"  =  26.  Using  Equation  (39),  E"  «  ~  2447.23cm“^. 

McClatchey  (Ref.  10:4)  provides  the  following  values 
of  Q^(T): 

Temp("K)  175  200  225  250  275  296 


Q..(T) 


1.0095  1.0192  1.0327  1.0502  1.0719  1.0931 


McClatchey  also  states  that  the  ratio  of  the  rotational 
partition  functions  varies  as: 

Vi>  ^  ( 

Qj.(T)  I 

where  j  =  1  for  CO2. 

This  expression  is  the  result  of  an  approximation  of  Q^. 
Instead  of  summing  over  all  values  of  j  in  Equation  (34), 
integration  results  in  the  approximation  (Ref.  38:125): 


Qr  = 


ahcB 


where  the  symmetry  factor  of  C02»  0*2.  Therefore, 
varies  as  the  temperature,  T.  The  U.S.  Standard  Atmosphere 
1976  (Ref.  39:53)  was  used  for  temperatures  and  pressures 
at  various  altitudes.  Using  the  values  for  Q^,  E”,  and 
Equations  (38)  and  (40),  "off "-resonance  cross  sections  for 
CO2  as  function  of  altitude  are  calculated.  Then,  to 
obtain  an  extinction  coefficient  as  a  function  of  altitude, 
it  is  necessary  to  multiply  each  cross  section  by  the 
molecular  number  density  of  CO2  at  the  corresponding 
altitude.  The  number  density  of  CO2  is  assumed  to  be 
uniformly  mixed  at  330pE«i  of  the  total  molecular  density  as 
a  function  of  altitude  (Ref.  23:6).  Absorption  cross 
sections,  number  densities,  and  extinction  coefficients  for 
"off "-resonance  CO2  can  be  found  In  Appendix  H. 


■  v:-.v 


Extinction  Coefficient  of  Carbon  Monoxide 


Calculation  of  extinction  coefficients  of  carbon 
monoxide  as  a  function  of  altitude  closely  parallels  that  of 
COj  above.  However,  a  few  differences  exist.  McClatchey 
(Ref.  10:4)  reports  that  does  not  vary  as  temperature. 
Therefore,  '^v  s  *  1.  Also,  in  Equation  (40),  j  «  1 

for  CO.  The  energy  of  the  lower  transition  state  (E”)  for  a 
diatomic  molecule  can  be  approximated  as  (Ref.  51:V-2) 

E5  „  =  B"J"(J"  +  1)  +  (V"  +  i)'Sr7  -X^(v"  +  i)^"Sr~  (42) 

J,  V  e  e  e 

where  for  CO  as  previously  stated 

B"  =  1.9225cm“^ 

=  2169. 74cm"^ 
e 

X  e  =  .0061 

The  lower  vibrational  energy  state  of  the  (1-0)  band  is 
ground  state,  so  v”  «  0.  The  rotational  line  closest  to 
"of f "-resonance  frequency  is  the  P(19)  line,  so  J"  »  19. 
Therefore,  E"  has  a  value  approximately  equal  to  1812.12cm'~^ 
From  Equation  (38),  absorption  cross  sections  as  a  function 
of  altitude  can  be  calculated.  As  previously  stated  CO  is 
considered  a  uniformly  mixed  gas  with  a  number  density  of 
.075ppm,  so  number  densities  at  various  altitudes  can  be 
calculated.  Appendix  I  lists  the  cross  sections,  number 
densities,  and  extinction  coefficients  for  CO  as  a  function 
of  altitude  for  the  "off "-resonance  frequency. 


These  calculations  are  similar  to  those  above. 

McClatchey  (Ref.  10:4)  reports  the  following  values  for 
as  a  function  of  temperature: 

Temp(K)  175  200  225  250  275  296 

Q^(T)  1.004  1.007  1.013  1.022  1.032  1.046 

The  ratio  of  the  rotation  partition  functions  is  given  by 
Equation  (40),  but  for  ozone,  j  «  1.5.  Since  the  spectra  of 
ozone  is  very  complicated,  with  numerous  lines  contributing 
to  absorption,  E"  is  taken  as  an  average  of  the  sum  of  the 
lower  energy  levels  of  the  170  evaluated  lines.  Flaud  and 
Camy-Peyret  (Ref.  37:192)  provide  values  of  E"  for  each 
Individual  absorption  line.  Calculating  the  mean  gives 
E”  »  827.400cm~^.  With  these  values,  and  Equation  (38), 
cross  sections  as  a  function  of  altitude  are  calculated. 

Ozone  is  not  a  uniformly  mixed  gas.  However,  U.S. 

Standard  Atmosphere,  1976  provides  number  densitites  as 
a  function  of  altitude  for  a  mid-latitude  ozone  model  (Ref.  39:38) 
Using  these  values  along  with  the  cross  sections,  extinction 
coefficients  as  a  function  of  altitude  can  be  calculated 
from  Equation  (17),  which  are  found  in  Appendix  J.  In 
addition,  the  extinction  coefficients  were  graphed  versus 
altitude,  and  interpolation  provided  the  additional  values 
found  in  Appendix  J. 


It  should  be  noted  that  the  extinction  coefficient 


of  ozone  does  not  vary  a  great  deal  with  altitude  due  to  a 
balancing  of  a  decreasing  cross  section  multiplied  by  an 
increasing  number  density,  which  is  in  direct  contrast  to  CO 
and  CO2  extinction  coefficients,  which  decrease  rapidly 
with  altitude.  Therefore,  a  check  to  determine  whether 
absorption  above  30  kilometers  is  significant,  where  Doppler 
broadening  is  dominant,  must  be  made.  From  Equation  (36) 
the  line  intensity  at  31  kilometers  of  the  line  contributing 
the  most  to  absorption  at  the  "off "-resonance  frequency 
(2057.024cm  is  given  by; 


S(T) 


s,T  ,rv!sL°Elvl 

““s’  Qy<T)  Q^(T)J 


exp 


1.439(T  -  T  )E" 
s 


TT. 


— 


where 


S(Tg)  =  5.77  X  10”^ ^cm”^cm^mol~^ (Appendix  G) 
T  s  227. 5 "K  at  31  kilometers. 

This  gives  S(T)  *  2.631  x  10”^^cm”^cm^mol'”^ 

Equation  (20)  gives  the  (Saussian  half-width  as 


^0/  \  * 

oijj  » - ((ln2)2kT/mj 

where  for  ozone 

m  ■  48amu  (1.6604  x  lO'^^^Vamu) 

'^o  ■  2057.054cm”^  (Appendix  G) 

_  -  1.60  x  lO'^cm"^. 


This  gives  oc 


The  cross  section  is  given  by  Equation  (20): 


cr(\))  *  (S/ac_)  (In2/Tr)*  exp  -(\)  -  \J '‘ln2/cx- 


Substituting  the  above  values  gives  2  0.  Therefore, 


extinction  above  30  kilometers  can  be  neglected. 


Extinction  Coefficient  of  Water  Vapor 

These  coefficients  are  calculated  in  the  same  manner 
as  the  ozone  above  with  only  minor  differences.  McClatchey 
(Ref.  10:4)  reports  that  the  vibrational  partition  function 
does  not  vary  with  temperature.  Since  H2O  has  the  Scune 
structure  as  ozone,  the  value  of  j  in  Equation  (40)  is  also 
1.5.  Averaging  over  the  107  lower  energy  states  provided  by 
Plaud  and  Camy-Peyret  (Ref.  36)  gives  a  mean  E"  of 
1709. 79cm""^.  Using  these  values  and  Equation  (38), 
absorption  cross  sections  as  a  function  of  altitude  are 
calculated. 

The  number  density  of  water  vapor  is  not  uniformly 
mixed.  Instead,  it  varies  a  great  deal,  depending  on  the 
atmospheric  conditions.  U.S.  Standard  Atmosphere,  1976 
provides  data  for  a  mixing  ratio  of  water  vapor  mass  to  dry 
air  for  an  average  mid-latitude  model  for  various  altitudes 
(Ref.  38:44).  After  converting  this  data  to  molecular 
number  density  by  volume,  extinction  coefficients  as  a 
function  of  altitude  can  be  computed  using  Equation  (17), 
results  are  found  in  Appendix  K. 


Calculation  of  "Of f" -Resonance  Return  Signal  Power 


The  returned  signal  power  for  the  "of f "-resonance 
frequency  is  determined  from  the  basic  lidar  equation, 
Equation  (4): 


/  \ 

1 

— 

f  TC  1 

-2  h  *  "cOj 

m  m 

23  R.2  I 

3 

where 

=  2  joules/lOOn  sec 
eK  =  ( .1)  ( .3)  =  .03 

=  r^  =  7.854  X  10“^ km^ 
=  R  =  1km 


To  determine  the  returned  power  from  various  altitudes  the 
exponential  of  Equation  (43)  must  be  evaluated.  Since  the 
extinction  coefficients  do  not  vary  as  simple  functions  of 
r,  numerical  integration  is  used  where: 


/o^^  ^2  **C02  ^  ^  ^2  **C02  ^2 


>Ari 


(4' 


where 


0,  H,0 

^2  ’  ^2  ^2  ^2R  ^2M 


Ar^  ■  1km 


where  ^2^  ^2M  *  extinction  coefficient  due  to 


Rayleigh  and  Hie  scattering. 


As  an  example,  consider  an  evaluation  of  the  return 


signal  from  a  one-kilometer  cell  from  four  to  five 

kilometers  in  altitude,  or  when  R.  =  196km.  For  this 

3 

example.  Equation  (44}  takes  the  form  of: 


196 


(  ?  +N 


170 


CO, 


,)dr 


196  H.O 

:  H  <^1* 

i*170  ^ 


CO 


2R 


2M 


+  N 


where  summations  are  done  over  extinction  coefficients  found 

in  the  applicable  appendices.  In  this  case, 

196 

f  (  52  +  Sco  -  •282*- 

170  ^ 

The  only  additional  unknown  term  in  Equation  (43)  is  the 
total  volumetric  backscattering  coefficients,  which  is  the 
sum  of  the  Rayleigh  and  Mie  backscattering  coefficients. 
Equation  (15)  gives  the  expression  for  the  Rayleigh 
backscattering  coefficient  as: 


ijR  8  IT  ^  iR 

The  direction  of  Mie  scattering  depends  heavily  on  particle 
size,  but  is  only  generally  considered  isotropic 
(Ref.  24:1502).  Therefore,  the  backscattering  coefficient 
due  to  Mie  scattering  is  given  by: 


Using  the  scattering  coefficients  from  Appendices  A  and  B  at 
196]an,  the  total  backscattering  coefficients  can  be  calcu¬ 


lated  from: 


2,  j 


2.32  X  10 


-4 


47r 


+  1-28  X  10  ^(3)  ^  ^  10"^km"^sr" 


87r 


Substituting  these  values  into  Equation  (43)/  the  returned 
signal  power  from  scattering  from  an  altitude  of  four  to 
five  kilometers  (R^  ®  196)  is: 

P2j  =  1.27  X  10"^°  Watts 


Similar  calculations  are  done  for  Rj  =  190  to  200. 

Results  of  the  total  backscattering  coefficients  and 
returned  signal  power  at  various  altitudes  are  displayed  in 
Table  I. 


Background  Noise  Calculations 

Three  basic  sources  of  background  exist: 

1.  Radiance  from  the  earth. 

2.  Radiance  from  the  sun. 

3.  Radiance  from  the  sky. 

Radiance  from  the  earth  is  considered  as  blackbody  radiance 
at  a  temperature  of  296 *’K.  The  equation  governing  the 
cunount  of  power  in  watts  received  at  the  lidar  detector  due 
to  the  earth's  background  radiance  is  (Ref.  20:25): 


Table  I 

Backscattering  Coefficient  and  Returned  Power 
at  "Of f "-Resonance  Frequency 


Rj  (km) 


fl2j(km- 

•^sr-1) 

2.520 

X 

10-s 

1.655 

X 

10-« 

2.867 

X 

10-« 

5.195 

X 

10-6 

7.516 

X 

10-6 

1.120 

X 

10-6 

1.861 

X 

10-6 

3.717 

X 

10-6 

3.082 

X 

10“^ 

5.135 

X 

1 

o 
1— 1 

7.960 

X 

10-* 

P2j(W) 


3.225 

X 

10 

2.068 

X 

10 

3.457 

X 

10 

5.918 

X 

10 

7.793 

X 

10 

9.935 

X 

10 

1.271 

X 

10 

1.580 

X 

10 

5.932 

X 

10 

2.876 

X 

10 

3.986 

X 

10 

where 


W(X)  «  Spectral  radiance  of  the  earth 

-2  -1  -1 
(Wcm  p  sr  ) 

AX  =  Spectral  band  width  (p) 

D  =  Diameter  of  receiver  telscope  (cm) 

©  *  Field  of  view  of  the  telescope  (assumed 

conical  with  0  full  cone  angle) 

K  =  Transmittance  of  receiver  optics 

T'  s  Atmospheric  transmittance 

For  operation  on  the  "off "-resonance  frequency,  T'  can 

Si 

be  determined  from: 


T '  =  exp 
a 


c-  c 


‘  ^2  *  "cOj 


] 


(47) 


Extinction  coefficients  are  found  in  Appendices  A,  B,  and  H-K 

-2 

In  this  case  T'  =»  6.52  x  10  .  Also,  K  is  again 

a 

taken  as  .3.  To  minimize  background,  ©  and  X  should  be 

minimized.  ©  is  selected  to  minimize  the  background 

signal  from  the  target,  while  at  the  same  time  ensuring  that 

the  FOV  encompasses  the  laser  spot  size  at  the  target.  © 

-4 

is  assumed  to  be  10  radians  (Ref .  20:28),  which  gives  a 

spot  size  on  the  earth  of  diameter,  d  »  ©R^  «  200  x 
3  -4 

10  (10  )  »  20m.  If  the  transmitted  laser  beam  is 

designed  to  also  have  a  spot  size  of  dieuneter,  d  «  20m,  then 
the  spot  size  (radius)  leaving  the  shuttle  is  given  for  a 
far-field  approximation  by  (Ref.  49:4-13): 


where 


w(z)  =  Spot  size  at  range  R^(m) 

Substituting  applicable  values  gives  w^~  3.09cm,  a 
reasonable  value  for  this  system. 

The  width  of  a  CO2  absorption  line,A\), 

-1  1  -4 

averages  .07cm  ,  or  since  X  =  ^/AX  =  1.65x10  ^Jl 

• 

For  this  study  a  spectral  bandwidth,  AX  =  6.75  x  10~^ji,  was 
selected  based  on  available  filters  at  about  5^.  The 
radiance  of  the  earth  is  given  by: 

Lg  =  W(X)  AX 

Since  the  earth  radiates  as  a  nearly  perfect  blackbody  at 
X  -  4-5fi  the  radiance  can  be  determined  from 
(Ref.  41:24): 


where 

H  V  Y  (emittance  in  spectral  window  of 
a  perfect  blackbody) 

*  /  M^(X)dK  -  /  M^(X)dX  (50) 

0  ®  0  ® 

€  (Emissivity)  ■  .85  (Ref.  29:3-90). 


'-•'v 


Radiation  tables  (Ref.  40)  provide  ratios  of  the  emittance 
of  a  blackbody.  The  total  emittance  is  given  by  the 
Stefan-Boltzmann  Law: 


M 


e 


aT 


(51) 


where 


a 


5.672  X  10"®Wm~^ 


To  get  aX  =  6.757  x  10  the  minimum  and  maximum 

wavelengths  used  were  =  4.858108p  and  X2 

=  4.864865ti.  Using  Equations  (49),  (50),  and  (51),  =  5.749 

-7  -2  -1 

X  10  Wcm  sr  .  Substituting  appropriate  values  into 
Equation  (46),  gives  P®  =  8.17  x  10~^^W. 

Calculation  of  background  due  to  the  sun's  radiation 
is  very  similar.  The  power  received  is  given  by: 


psun 

B 


W(X)  X  3;|-  9  T  ^r  K 


(52) 


where 

r  3;  reflectivity  of  earth  at  4.86p  estimated  at  .15 
(Ref.  29:3-90) 

Atmospheric  transmittance  is  squared  since  radiation  must 

pass  through  the  atmosphere  twice.  The  sun  is  also 

considered  a  blackbody,  but  at  a  temperature  of  5900”K. 

The  radiance,  L^,  is  calculated  using  radiation  tables  in 

a  similar  manner  to  that  described  above.  L  is 

e 


t  »  VV.  ^2 

calculated  as  x  >er  Wcm~'^sr  .  Prom  Equation  (52), 

the  background  noise  due  to  the  reflected  sun's  radiance  is 

3.^5-  ic'L'' 

calculated  as  x  J^^W. 

0 

Background  noise  as  a  result  of  sky  radiance  is  more 
difficult  to  analyze.  This  background  is  a  result  of  scat¬ 
tering  of  the  sun's  radiance  and  emission  by  atmospheric 
constituents.  However,  essentially  the  scattering  occurs 
below  34  and  the  thermal  emission  above  44  (Ref.  29:3-71). 

To  account  for  this  additional  background  the  experimental 

data  by  Bell  (Ref.  42:1318)  has  been  used.  Bell  estimates 

-4  -2  -1  - 

the  sky  spectral  radiance  at  4.864  as  2.13  x  10  Wcm  4  sr 

—6  —2 

Multiplying  by  aX  gives  a  radiance,  =  1.439  x  10  Wcm  sr 
Losses  due  to  additional  atmospheric  scattering  are  ignored 
so  the  background  due  to  sky  radiance  received  at  the 
detector  is  given  by: 


=  2.66  X  10"^^W. 


(53) 


The  total  background  is  the  sum  of  the  three  individual 
background  noises: 


P  =  P®  +  P®^"  +  P®’'^ 
*  3^  X  i^w 


(54) 


P 


Siqnal-to-Noise  Calculations  (Direct  Detection) 

The  detector  offering  the  most  efficient  normalized 
detectivity  in  the  4.86pt  region  is  the  indium  antimonide 
(InSb)  photovoltaic  infrared  detector  operating  at 
Kildal  and  Byer  state  that  a  photovoltaic  InSb  detector  is 
always  dark-current  limited  for  normal  bandwidths.  In  this 
case  they  give  the  normalized  detectivity  as  (Ref.  43:1648): 


V2e<VV 


or  by  rewriting,  the  dark  current  is  given  as: 

1 

^  2eD*^  (56) 

where  for  a  Ford  Aerospace  and  Communications  InSb  infrared 
detecter  (ISC-386  series) 

JV 

Sjj  (Sensitivity  of  the  detector  =  2.0  /W 
e  =  1.602  X  10“^^C 

D*  »  8.6  X  10^°cm  Hz^w"^  for  a  0^  »  60* 

However,  the  field  of  view  developed  earlier  has  9^  ■ 

5  X  lo'^rad.  Therefore,  it  is  necessary  to  develop  a  new 
normalized  detectivity  from  the  equation  below  (Ref.  44:45): 


D**  a  D*sin0| 

D**  »  7.45  X  10^°cmHz^sr*w"^. 

Therefore,  for  a  FOV  of  10”^rad,  D*  =*  1.49  x  lO^^cmEz^w”^. 
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Substituting  into  Equation  (56)  i^  *  1.59  x  10  A. 


However,  the  current  due  to  the  returned  laser  signal  is 
approximately : 


i„  =»  (1.27  X  10“^°)  (2.0)  at  R.  =  1961cm 
s  D 

i  =  2.54  X  10“^°A 
s 

Therefore,  since  the  FOV  is  so  narrow,  the  InSb  detector 
cannot  be  considered  dark  current  limited,  and  i.  can  be 

U 

neglected  in  future  calculations. 

The  power  signal-to-noise  ratio  for  photovoltaic 
photodiode  is  given  by  (Ref.  41:169): 


i|  =  Detector  quantum  efficiency 

f  =  frequency  of  operation  corresponding  to  "on" 
or  "off"  resonance  (Hz) 

Af  =  electrical  bandwidth  of  the  detection  system  (Hz) 
m  s  Order  of  the  modulated  optical  beam. 

=  Effective  temperature  of  the  load  resistor, 
including  cunplifier  noise. 

=  Load  resistor  of  the  detection  system  (Rj^  «  R^) 

For  an  "of f "-resonance  pulse,  f  =  6.1668  x  lO^^Hz, 
r|  is  talcen  as  a  conservative  .3  (Ref.  20:264),  and  m  »  1. 

The  electrical  bandwidth  must  be  large  enough  that  it  is 
responsive  to  the  return  laser  signal.  Commonly,  the 
minimum  bandwidth  is  taken  as: 


where 

T ( intergration  time  of  the  returned  power  signal)  » 

^  *  6.67  X  10"®sec 
c 

4 

In  this  case,  Af  =  7.495  x  10  Hz. 

The  effective  temperature  (T^)  includes  the 
effective  noise  temperature  of  the  eunplifier,  T^.  T^  is 
determined  from  the  2uaplifier  noise  figure,  ?,  where 
(Ref.  41:169):  [ _ 


If  P  has  a  realistic  value  of  2  (Ref.  41:170),  then 
T-  *  290 “K.  The  effective  temperature,  T  ,  is  the  sum 

A  6 

of  the  detection  system  temperature  plus  the  amplifier 

temperature : 

*  T  + 

If  T  =  290 “K,  then  =  580 “K. 

Topically,  for  a  photodiode  the  load  resistor  is  selected  to 

be  less  than  the  impedance  of  the  detector,  R^,  which  for 

ISC-386  series  photodiode  is  25  x  10^  ohms,  and  much 

greater  than  the  spread  resistance,  R  ,  typically  50  ohms. 

s 

For  this  study,  R^  ^  1000  ohms. 

As  an  example,  the  signal-to-noise  ratio  for  a 
return  signal  from  R^  =  196)an  is  given  by: 


C 

As  can  be  seen,  (  /1H) 2  does  not  meet  the  minimum  require- 
S  2 

ment  of  (  /N)  »  2  as  given  in  Equation  (12).  Clearly, 

the  dominating  noise  (second  term  in  the  denominator)  is  the 
Johnson  or  thermal  noise  of  the  detector,  load  resistor,  and 
amplifier.  Unfortunately,  there  is  little  that  can  be  done 
to  reduce  the  thermal  noise  the  two  orders  of  magnitude 

C 

necessary  to  achieve  an  acceptable  (  /N)2> 

Heterodyne  Detection 

The  best  alternative  to  achieve  an  acceptable 

C 

(  /N)2  without  changing  any  of  the  given  parameters  of 

63 


the  study  is  by  coherent  or  heterdyne  detection,  where  the 
incoming  wave's  amplitude  and  phase  are  converted  to  a  low 
frequency  in  the  rf  or  microwave  region.  Heterodyne  detec¬ 
tion  was  initially  proposed  for  IR  or  visible  frequency 
signals  in  the  19€0's,  and  has  since  successfully 
demonstrated  an  numerous  occasions  including  recent  experi¬ 
ments  with  pulsed  CO2  lidars  by  Cruickshank  (Ref.  45:290) 
and  Lundquist  (Ref.  46:2534). 

To  obtain  this  low  frequency  signal,  the  incoming 
signal  is  mixed  with  that  of  a  continuous  beam  of  slightly 
different  frequency  from  a  local  oscillator  (See  Figure  4). 
The  effects  of  this  mixing  are  discussed  below  for  the 
simple  case  of  equal  phase,  constant  aunplitude  co-planar 
waves  (Ref.  47:24).  The  total  electric  field,  E,  of  the 
incident  waves  is  given  by: 

E  =  Ej^qCos  +  Ej^jCOs  Wj^jt  (61) 

where 

E^q  *  Electric  field  of  the  local  oscillator  with 
frequency, 

Ej^j  «  Electric  field  of  the  returned  signal  at 
frequency, 

The  photodector  current  from  the  combined  signal  is  given  by: 
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Figure  4; 


Block  Diagram  of  a  Laser  Radar 
With  Heterodyne  System  (Ref.  48:258) 
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The  power  can  also  be  expressed  as: 


E^(t)A. 


P{t) 

where 


(63) 


=  Detector  area 
-  Impedance  of  free  space 


Combining  Equations  (61),  (62),  and  (63): 

^  I*  ''to  "  “to‘>  +  ‘  “’2  “lj‘' 

*  -  “to’'  *  “to’'] 

The  final  term  can  be  ignored  since  the  elctronics  will  not 

resolve  the  combined  frequencies.  Recalling  the  i^^^  ® 

en  A .  M  0 

ri-CT  e:„  and  i^^  *  Ef ^  gives: 


2z  hf  "LO 
o 


i(t)  +  ij^Q  +  iij  +  2  rwlj  ‘  ^  “ij  "  “lo’^ 

C 

To  calculate  the  power  (  /N)^^  due  to  the  intermediate 
frequency  (i.f.)  current,  take  the  mean-square  i.f.  current: 

^if  “  ^  ^LO^Tj^  *  ^^LO^ij  (66) 


If  the  shot  noise  is  dominated  by  the  local  oscillator  beaun, 
by  design,  then  the  mean-square  shot  noise  current  for  a 
photodiode  is  given  by  (Ref.  47:25): 


where 


=  Bandwidth  of  the  i.f.  channel  following  the 
detector  (Hz) 

Including  the  thermal  noise  of  the  detection  system  gives 
the  power  signal-to-noise  ratio  as  (Ref.  47:41); 


Combining  Equations  (66),  (67),  and  (68)  gives 


r 

^ij 

hfB^j  1 

f  2  ^  2kT^hf  1 

i  «r.  1 

(69) 


In  order  to  ensure  that  the  shot  noise  of  the  local 
oscillator  dominates  the  thermal  noise  of  the  detection 
system: 


^LO  >> 


21cT  hf 
e 


(70) 


Assuming  the  same  values  for  the  variables  as  previously 
given,  the  local  oscillator  power  must  be  greater  than  4.249 
X  10~^W.  A  reasonable  value  would  be  .IW,  accounting  for 
the  required  frequency-doubling  crystal.  If  a  local 
oscillator  is  used  then  the  power  signal-to-noise  ratio  is 
given  by: 


Several  important  features  can  be  summarized  at  this  time 
(Ref.  47:26): 


1.  The  i.f.  signal  power  is  proportional  to  the 
returned  signal, 

/  >  ^  j 


inversely  proportional  to 


the  bandwidth. 

3.  The  output  i.f.  spectrum  for  a  broad-band  signal 
is  an  exact  replica  of  the  input  optical 
spectral  distribution,  provided  the  local 
oscillator  is  a  single-fequency  wave  with 


nonfluctuating  phase. 
Table  II  lists  values  for/ 


for  values  of 


P_.  with  values  of  j  =  190  to  200Jcm,  assuming  B.  .  *  Af. 

/  e\ 

Equation  (71)  gives  the  common  equation  for  (  ~j  .  However, 
two  problems  exist  with  this  analysis.  First,  it  does  not 


account  for  random  signal  flucatuation,  and  second,  the 
coherence  of  the  signal  may  not  be  maintained  during  its 
analysis.  In  a  pulsed  signal,  the  coherence,  at  best,  is 
usually  the  pulse  length.  If  Bj.£t'<l,  then  a  temporal 
averaging  factor  must  be  added.  A  more  accurate  expression  of 
has  been  developed  by  Elbaum  (Ref.  48:257): 

/s)  . 

I  Tvar 


where 


and^i(t)^'^=  average  values  of  output 
current  in  the  presence  and  absence  of 
signal,  respectively. 

For  a  photodiode,  this  takes  the  form  of  (Ref.  22:1175): 


hi' 


(s\ 

I** /Hi  ^  ^^^1  +  j  t 


k(T+T,) 

a 


where 


ProM®\^  1 

G(heterodyne  conversion  gain)  =  — -I  | - s - ; 

2G_  V  hf  /  1  +  (  /fe)' 


where 

Gjj  *  Incremental  shunt  conductance, 

f^  =  Rolloff  frequency  of  the  photodiode. 

If  P--  shot  noise  again  is  designed  to  dominate 

liU 

other  noise  forms,  then  Equation  (73)  becomes: 


<8H>]  * 


hfB^  *  ‘■ij 


Note,  that  the  term  1 (^  /c) (B 


must  be  at  least  four 


if  is  the  dominate  term  in  the  denominator,  an< 


is  to  be  greater  than  two.  Therefore,  select  »  10  Hz. 
Table  III  lists  values  of|§)  for  R.  «  190-200km 

,  r/H2  J 


I  for  R. 
fH2  J 


190-200km 


with  B, 


10  Hz.  Note  that  all  values  exceed 


the  minimumf^l  requirement  of  two  given  in  Equation  (12). 


Detection  of  Hot  C0-»  in  the  Atmosphere 


It  is  conceivable  that  shuttle  observations  could 
detect  CO2  hot  spots  in  the  atmosphere  due  to  a  high 
energy  laser  facility  or  fossil  fuel-burning  plant. 

Equation  (4)  gives  the  basic  DIAL  equation  for  the  return 
power  from  a  shuttle  lidar.  If  a  quantity  of  hot  CO2  gas 
is  introduced  into  the  atmosphere  then  the  basic  equation  is 
altered  to: 


"ii  = 


-  fl.  .e*pi  -2  (  q  + 


a,  +  a ,^)dr 


where 


=  Number  density  of  hot  CO2  molecules. 

=  Absorption  cross  section  for  CO2  at  a  given 
temperature,  T. 

This  analysis  will  determine  the  maximum  number  density  of 
11 

hot  C02/  N|2q  /  that  can  be  detected,  and  still  meet  the 
minimum  ^  requirement  of  two  as  stated  in 


requirement  of  two  as  stated  in 


Equation  (12).  The  minimum  "off "-resonance  power  is 
determined  from  Equation  (74): 


R)* 


2hfB, 


2j  min 


4.42  X  10 


In  order  to  determine  the  number  density  of  the  hot 


CO2  that  will  give  that  minimum  return  power,  it  is 
necessary  to  determine  the  absorption  cross  sections  at 
various  high  temperatures.  Equation  (38)  gives  the 
"of f "-resonance  absorption  for  CO2  as  function  of  pressure 
and  temperature: 


(P,T) 


a  ■ 

2CO2 


Qv^^s> 

Q^{T) 


VV 

Qr(T) ^ 


exp 


1.439(T-T  )E 
_ s 

TT_ 


P_ 

P. 


where  <^200  “  1*361  x  lO'^^cm^,  and  E"  =  2447.23cm~^. 

Unfortunately,  McClatchey  does  not  provide  values  for  the 
vibrational  partition  function,  Q^,  for  temperatures  over 
296*K.  Therefore,  it  is  necessary  to  compute  values.  As 
previously  stated  the  vibrational  partition  function  is 
given  by; 

E. 


Qv  =  2-  exp 

Vi,V2,V3,l 


[- 


A  good  approximation  is  the  harmonic  oscillator 
approximation,  which  assumes  the  energy  levels  are  integer 
multiples  of  the  fundamental  frequencies.  This,  along  with 
the  identity  ^  exp  *  (1-exp^E 


gives  the  harmonic  approximation  of  the  vibrational  partition 
function  (Ref.  38:52): 


Qv  - 


(  ^  ^ 

II 

2 

f  1  \ 

1 

^  1 

(77) 


where 

E^/k  *  1.439Vj^ 
E^/k  *  1.439^2 
E^/k  =  I.439V3 


=  1330.0cm' 
"^2  “  •  3cin 

~  2349.3cm' 


Therefore,  the  ratio  of  the  vibrational  partition  function 
is  given  by: 


Q^(T) 


-1.439\) 
(1-e _ 

-1.439^0 

(1-e 


,  -1.439\)  ,  ,  -1.439\)  , 

^/T)(l-e _  /T)^  (1-e _ VT) 

,  -1.439Sr  o  ,  -1.439\), 

/r)  (1-e  /'SJ  (1-e  ^/T) 

s  s  s 


The  ratio  of  the  rotational  partition  functions  is 
again  given  by  Equation  (40): 


k!Z2> 

^(T)  “It  j 


where  j  =  1  for  CO2. 

Four  assumptions  are  made  in  this  analysis: 

1.  The  additional  hot  CO2  does  not.  significantly 

add  to  the  atmospheric  pressure  as  a  function  of 
altitude. 


2.  The  gas  has  diffused  and  cooled  above  2 
kilometers . 


3.  The  laser  illuminates  a  uniform  area  of  hot  CO 


2* 


4. 


The  hot  CO2 
cell ,  and  N 


occupies  only  a  single  one-kilometer 


__  ,  and  are  constant  over 


the  cell. 

With  these  assumptions,  and  pressure  data  from  the 
U.S.  Standard  Atmosphere,  1976,  "off "-resonance  absorp¬ 
tion  cross  sections  can  be  calculated  for  temperatures  from 
T  =  400“K-1000“K.  Appendix  L  lists  these  values  for 
altitudes  from  0-2  kilometers. 

By  combining  the  basic  DIAL  equation.  Equation  (4), 
and  the  modified  DIAL  equation  that  accounts  for  the  hot 
CO2,  Equation  (76),  relationship  between  the  returned 
signals  can  be  expressed  as: 


(79) 


The  minimum  acceptable  returned  power,  P^.  ,  is 

substituted  for  to  get  the  maximum  .  Values  for 
P2j  are  found  in  Table  II.  Therefore,  for  R2j  =  198-200km, 
the  intsrgrated  optical  thickness  of  hot  CO>  is  given  as: 


By  using  a  one  kilometer  cell,  and  assuming  constant  density 
cross  section  over  the  cell,  as  stated  in  assumptions  (3) 
and  (4),  the  maximum  hot  CO2  density  can  be  computed. 

Results  as  a  function  of  temperature  for  =  198-200  are 
found  in  Table  IV. 


Table  IV 

Maximum  Hot  C0„  Concentration  fo 


Temp  (®K) 


198km 


199km 


200km 


H  I 

^  2CO2 

Ccm^> 

1 - 

1.334 

X 

10-20 

4.824 

X 

10-20 

1.017 

X 

1.581 

X 

10-^^ 

2.055 

X 

10-15 

2.382 

X 

10-15 

2.546 

X 

10-15 

1.507 

X 

10-20 

5.541 

X 

10-15 

i 

1.149 

X 

10-15 

1  1.787 

X 

10-15 

i  2.322 

X 

10-15 

2.692 

X 

10-15 

2.876 

X 

10-15 

1.699 

X 

10-20 

6.145 

X 

10-20 

1.295 

X 

10-15 

2.015 

X 

10-15 

2.617 

X 

10-15 

3.035 

X 

10-15 

3.243 

X 

10-15 

2.70  X  10 
7.46  X  10 
3.54  X  10 
2.28  X  10 
1.75  X  10 
1.51  X  10 
1.41  X  10 


2.14 

X 

10 

5.91 

X 

10 

2.18 

X 

10 

1.80 

X 

10 

1.38 

X 

10 

1.20 

X 

10 

1.12 

X 

10 

1.32 

X 

10 

3.66 

X 

10 

1.74 

X 

10 

1.12 

X 

10 

8.60 

X 

10 

7.42 

X 

10 

6.94 

X 

10 

Error  Analysis  of  DIAL  Calculations 


Equation  (4)  gives  the  returned  signal  power  for  a 


range-resolved  DIAL  system; 


eP^KA^  /  ct\  r  R.  ■) 

>ij  <ei  *  •'coj  "i  >*-^3 


When  i  =  1  designates  the  "on "-resonance  frequency,  i  =  2 
designates  the  "of f "-resonance  frequency,  j  =  1  designates  a 
range  cell,  and  j  =  2  designates  the  j  +  1  range  cell,  four 
equations  of  the  basic  DIAL  equation  can  be  expressed. 
Manipulation  of  these  equations  gives  the  expression  for  the 
optical  thickness  for  a  range  cell  to  R2: 


^11^22 

Aa  dr  »  i  In 

\  2  P,2P21 


®12®21  ^2 

In  -2  f  (  C  ,  -  C 

«11«22  «1 


where 


Ac  _  a 


Since  the  volumetric  backscattering  coefficient  varies 
slowly  in  infrared  as  a  function  of  frequency,  it  can  be 
assumed  that  8^2  -  ^21  ^1  -  ^2*  Inhere  fore, 

the  second  term  of  Equation  (80)  can  be  neglected.  Also, 
Korb  in  his  analysis  selected  his  "on"  and  "off "-frequencies 
based  on  the  analysis  that  extinction  coefficients  for  other 
than  CO2  were  equal  for  the  two  frequencies.  For  that 
reason,  the  third  term  of  Equation  (80)  can  be  neglected. 


Finally,  if  it  is  assumed  the  number  density  and  cross 
section  are  constant  over  the  range  cell,  then  the  number 
density  can  be  expressed  as: 

Range-resolved  differential  absorption  lidar  data 

reduction  involves  the  comparison  of  large  amplitudes  having 

small  differences.  The  error  of  the  CO2  concentration 

depends  heavily  on  errors  in  determining  the  magnitude  of 

the  returned  signals.  Errors  can  arise  from  laser 

wavelength,  background  noise,  aerosol  concentration, 

accuracy  of  H2O  vapor  profile,  lidar  calibration 

procedures,  and  transmissivity  due  to  atmospheric 

constituents.  Because  of  these  uncertainties  it  is  not 

adequate  to  describe  DIAL  measurement  accuracy  based  soley 

on|^  I  ,  but  an  expression  must  be  developed  that 
V”  /  Hi 

includes  the  random  coherent  noise  within  the  signal 
(Ref.  21:557). 

The  error  propagation  by  least  squares  ( variance ) 
for  a  function  f(x^.  .  .x^)  is  given  by  the  expression 


(Ref.  19:8): 


where 


2  ^ 


-•  f  ( X  X  ^ 

o  01  on 


where 


If  It  IS  assumed  that  the  absorption  cross  section  can  be 
accurately  verified  through  laboratory  experimentation,  and 
that  the  range  cell  height  can  be  accurately  determined  from 
data  analysis,  then  the  standard  deviation  reduces  to: 

,  r/  ,  \2  2,  2_  n 


i™)'  k  k  lif-] 


But  cannot  be  directly  measured.  Instead  P^^  can  be 


expressed  as: 


where 


Q  ®  Amplification  factor  (  /W) , 


^i(t)^  and^i(t)^^  =  Average  values  of  the  output 

current  in  the  presence  and  absence  of  signal  as 
stated  in  Equation  (72). 

Therefore,  the  variance  of  can  be  expressed  as; 

Combining  with  Equations  (84)  and  (85)  gives: 


Hi  j^Var  { 

Combining  with  Equation  (87)  gives 


ml 


are  measured,  then  Equation  (89)  has  the  final  form  of 
(Ref.  21:558); 


The  number  of  pulse-pairs  is  related  to  the  desired 
horizontal  resolution.  For  CO2  studies  the  NASA  Working 
Group  suggested  x  =  100km  (Ref.  7:18).  If  the  shuttle 
velocity  is  8km/ sec  (Ref.  15:626)  with  a  laser  repetition 
rate  of  lOHz  (Ref.  22:1181),  then  n  »  125  shots. 


Therefore,  Equation  (90)  can  be  expressed  as: 


But  the  "on"  and  "of f "-resonance  powers  are  related  by 
Equation  (5): 


(92) 


(93) 


Combining  Equations  (91),  (92),  and  (93)  gives: 
<^N  = 


_ ^\2  1 

/l2  AoA  R 1  n 


r'  1  R, 

—  (1  +  exp  4  /  N_^Aadr) 


S\2 
H22 


It): 


^  CO. 
O  2 


+  ,  3^2  (1  +  exp  4  1I(,Q  4odr) 

(5)h21  “  " 


(94) 


Prom  Appendices  C  and  D,  it  can  be  determined  that: 


14 . 91 CJ2  = 

Therefore,  Aa  =*  13.91  02*  Substituting  this  result  into 

Equation  (94),  and  using  extinction  coefficients  for  "off"- 

resonance  CO2  in  Appendix  H,  and^^j^  values  in 

Table  II,  the  standard  deviation  of  the  concentration  in 
-3 

cm  and  ppm  can  be  calculated.  These  results  are  given 
in  Table  V.  Clearly,  this  shows  that  even  above  1951an  that 
the  uncertainty  of  the  measured  concentration  may  preclude 
shuttle-borne  lidar  operation  at  the  proposed  frequencies, 
even  though  adequate  signal -to-noise  ratios  are  available. 


These  results  do  not  generally  compare  to  Korb  (Ref.  17) 
From  the  spectrum  in  that  abstract,  the  optical  depth  for 
the  "on"-resonance  frequency  for  a  two  kilometer  horizontal 
path  at  sea  level  is  given  as: 

?  5 

K  =  1  =  J  ^  1  *  10 

o  2 

For  a  sea  level  concentration  of  CO2  of  8.405 
15  -3 

X  10  cm  ,  the  "on "-resonance  cross  section  is  solved 

-22  2 

as  a  *  5 . 95  X  10  cm  .  This  does  not  agree  well 

with  the  previously  calculated  value  of  =  2.029  x 
-20  2 

10  cm  .  If  Kerb's  value  of  is  correct,  then  the 
error  of  measurement  will  be  much  less. 

The  fact  that  the  error  of  the  CO2  concentration 
is  so  large  is  due  to  the  large  optical  depth  for  the  "on"- 
resonance  frequency.  This  agrees  with  Remsberg  and  Gordley 
(Ref.  15),  who  predict  that  for  an  optical  depth  of  greater 
than  four,  the  on-line  signal  from  the  range  cell  will  be 
small  and  difficult  to  measure. 


Table  V 

Standard  Deviation  of  at  Various  Altitudes 


R2(kn) 

Aa  (cm^) 

-3) 

oN(  cm 

191 

-22 

1.67  X  10 

8.596  X  10^"^ 

-22 

14 

192 

2.87  X  10  ^ 

6.147  X  10^^ 

193 

-22 

4.83  X  10 

5.802  X  10^^ 

194 

-22 

7.94  X  10 

9.910  X  10^^ 

195 

-21 

1.27  X  10  ^ 

4.395  X  10^^ 

196 

-21 

1.99  X  10 

8.753  X  10^^ 

197 

3.07  X  10"^^ 

19 

3.391  X  10  ^ 

198 

4.63  X  10"^^ 

8.050  X  10^^ 

199 

-21 

6.87  X  10  ^ 

1.123  X  10^^ 

200 

-20 

1.00  X  10 

44 

5.261  X  10^ 

5.72  X  10 
1.99  X  10 
3.85  X  10 
5.36  X  10 

2.28  X  10 


VI.  Conclusions  and  Recommendations 


Conclusions 


On  the  basis  of  the  above  analysis  the  following 
conclusions  are  made: 

1.  A  range-resolved  differential  absorption  lidar 
is  the  only  system  capable  of  measuring  concen¬ 
trations  of  constituents  of  the  troposphere. 

Of  the  major  interfering  species  in  the  atmos¬ 
phere  at  4.86p  (carbon  monoxide,  water,  ozone, 
aerosols,  and  cardon  dioxide),  carbon  dioxide 
provides  the  major  contribution  to  the  total 
extinction  coefficient. 


2. 


3.  The  minimum 


(4 


required  is  determined 


to  be  two.  Direct  detection  could  not  produce 
this  result  due  to  the  thermal  noise  of  the  detec 
tion  system.  Heterodyne  detection,  however,  pro¬ 


duces 


(«): 


I  8  after  increasing  the 


bandwidth  to  10  Hz. 


4.  Addition  of  hot  CO2  below  2km  can  be  detected 

up  to  a  maximum  concentration  of  approximately 
15  -3 

10  cm  depending  on  altitude  and  temperature. 

5.  Error  analysis  of  measured  CO2  concentrations 
indicates  very  large  standard  deviations  particu¬ 
larly  below  five  kilometers.  This  is  due  to  the 
increased  extinction  of  the  laser  beam  due  to 


2/2 
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"on "-resonance  absorption.  This  does  not  agree 
well  with  Korb's  analysis  whose  "on "-resonance 
cross  section  is  computed  as  approximately  34 
times  smaller  than  the  cross  section  computed 
in  this  thesis. 

Recommendations 

Based  on  this  thesis  the  following  recommendations 
are  made: 

1.  Update  of  the  spectroscopic  data  of  the  4.82^ 
absorption  band  of  CO2  should  be  made  to 
determine  the  actual  line  intensities  and  widths 
so  exact  cross  sections  are  readily  available. 
This  would  determine  whether  the  P(34)  line  is 
acceptable  for  "on "-resonance  absorption,  or 
whether  another  absorption  line  or  band  in  com¬ 
bination  with  suitable  laser  would  provide  more 
meaningful  results. 

2.  Even  though  this  particular  investigation  did 
not  provide  suitable  results  to  recommend  a 
working  system,  it  appears  that  a  similar 
approach  at  another  frequency  in  the  4.82p 
spectral  region  could  track  CO2  concentrations 
with  relatively  small  error.  It  appears  the 
system  could  be  quite  effective  in  locating 
sudden  changes  in  CO2  concentrations  at  a 
large  industrial  complex  or  laser  facility. 
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Appendix  A 


Raleiqh  Scattering  Coefficients 
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Mie  Scattering  and  Extinction  Coefficients  (4.86u) 
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^2 

Absorption  Cross 

Section  ( "On "-Resonance ) 

ID 

-1  -1  -2 ) 

S  (cm  mol  cm 
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^ICOj 

(cm^) 

34) 
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lo" 
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C0~  Absorption  Cross  Section  ( "Off "-Resonance) 


V  (cm  )  S  (cm  mol 
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Total  a 
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CO  Absorption  Cross  Section  ( "Of f "-Resonance) 


Rot.  ID 


M.(cin"^) 
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HjO  Absorption  Cross  Section  ( "Off "-Resonance) 


\)_(cm"^) 
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10-24 

9.890 

X 

10 

2050.986 

.294 

X 

10-24 

8.214 

X 

10 

2050.649 

.106 

X 

10-23 

2.657 

X 

10 

2064.853 

.499 

X 

10-2^ 

8.923 

X 

10 

2064.854 

.150 

X 

10-20 

2.492 

X 

10 

2065.021 

.199 

X 

10-21 

3.170 

X 

10 

2065.644 

.156 

X 

10-24 

2.139 

X 

10 

2065.846 

.595 

X 

10-21 

7.787 

X 

10 

2049.730 

.723 

X 

10-24 

1.384 

X 

10 

2048.649 

.270 

X 

10-23 

3.921 

X 

10 

2047.725 

.201 

X 

10-24 

2.368 

X 

10 

2046.796 

.369 

X 

10-23 

3.593 

X 

10 

2046.516 

.126 

X 

10-21 

1.162 

X 

10 

2046.486 

.114 

X 

10"24 

1.046 

X 

10 

2066.614 

.349 

X 

10-24 

3.865 

X 

10 

2067.778 

.250 

X 

10-23 

2.202 

X 

10 

2068.896 

.159 

X 

2070.683 

.113 

X 

2043.951 

.401 

X 

2043.943 

.652 

X 

2043.698 

.142 

X 

2043.165 

.148 

X 

2041.513 

.396 

X 

2041.497 

.520 

X 

2041.289 

.156 

X 

2071.950 

.847 

X 

2072.003 

.403 

X 

2072.488 

.124 

X 

2072.541 

.167 

X 

2073.379 

.186 

X 

2073.990 

.526 

X 

2074.253 

.181 

X 

2074.689 

.145 

X 

2040.306 

.171 

X 

2040.016 

.859 

X 

2039.946 

.347 

X 

2039.292 

.145 

X 

2038.094 

.989 

X 

2037.507 

.113 

X 

2036.620 

.252 

X 

2057.629 

.106 

X 

10"'‘ 

1.149 

X 

10 

10-2^ 

5.356 

X 

10-30 

10-21 

2.390 

X 

10-2® 

10-2" 

3.881 

X 

10-29 

10-2" 

8.145 

X 

10-30 

10-2" 

8.489 

X 

10-30 

10-22 

1.677 

X 

10-27 

10-21 

2.197 

X 

10-2® 

10-20 

6.418 

X 

10-2® 

10-2" 

3.873 

X 

10-29 

10-2" 

1.830 

X 

10-29 

10-2" 

5.28 

X 

10-30 

10-22 

7.065 

X 

10-28 

10-23 

7.083 

X 

10-29 

10-23 

1.861 

X 

10-28 

o 

1 

Is) 

6.379 

X 

10-27 

10-2" 

4.733 

X 

10-30 

10-2" 

6.232 

X 

10-30 

10-2" 

3.025 

X 

10-29 

10-23 

1.212 

X 

10-28 

10-2" 

4.697 

X 

10-30 

10-2" 

2.811 

X 

10-29 

10-21 

3.022 

X 

10-27 

10-2" 

6.166 

X 

10-30 

10-2" 

3.119 

X 

10-30 

,«-24 

,«-29 

2076.030 
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24 


1.297  X  10 


2076.698 

.315 

X 

10 

8.289 

X 

10 

2076.973 

.483 

X 

10“^^ 

1.236 

X 

10 

2077.041 

.114 

X 

10-24 

2.898 

X 

10 

2053.530 

.194 

X 

10-24 

4.227 

X 

10 

2035.309 

.582 

X 

10-24 

1.257 

X 

10 

2035.149 

.785 

X 

10-24 

1.671 

X 

10 

2034.360 

.898 

X 

10-24 

1.781 

X 

10 

2034.191 

.134 

X 

10-24 

2.618 

X 

10 

2034.141 

.267 

X 

10-22 

5.194 

X 

10 

2034.054 

.445 

X 

10-22 

8.591 

X 

10 

2034.034 

.314 

X 

10-22 

6.051 

X 

10 

2032.290 

.123 

X 

10-22 

2.048 

X 

10 

2078.570 

.103 

X 

10-22 

2.260 

X 

10 

2079.932 

.376 

X 

10-22 

7.298 

X 

10 

2081.874 

.457 

X 

10*22 

7.538 

X 

10 

2085.494 

.141 

X 

10*22 

1.722 

X 

10 

2087.408 

.205 

X 

10*22 

2.262 

X 

10 

2027.025 

.67 

X 

10*22 

7.017 

X 

10 

2026.603 

.179 

X 

10-22 

1.970 

X 

10 

2023.030 

.421 

X 

10-22 

3.711 

X 

10 

2019.071 

.577 

X 

10-22 

4.080 

X 

10 

2018.338 

.117 

X 

10-20 

7.963 

X 

10 

2016.836 

.367 

X 

10-20 

2.315 

X 

10 

2016.799 

.122 

X 

10-20 

7.680 

X 

10 

2009.333 

.197 

X 

10-22 

8.823 

X 

10 

2007.701 

.407 

X 

10-22 

1.702 

X 

10 

2089.742 

.724 

X 

10 

6.889 

X 

10 

2090.023 

.676 

X 

10-22 

6.323 

X 

10-“ 

2090.102 

.623 

X 

10-2^ 

5.800 

X 

10-“ 

2090.107 

.208 

X 

10-21 

1.936 

X 

10-“ 

2097.368 

.478 

X 

10-22 

2.991 

X 

10-“ 

2100.433 

.543 

X 

10-22 

2.935 

X 

10-“ 

2106.347 

.210 

X 

10-22 

8.793 

X 

10-“ 

2107.547 

.210 

X 

-22 

10 

8.380 

X 

10-“ 

2114.427 

.614 

X 

10-“ 

1.898 

X 

10-“ 

2115.018 

.234 

X 

io-“ 

7.087 

X 

10-“ 

Total 

° 2E^O  * 
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Appendix  G 


0, Absorption  Cross  Section  ( "Off "-Resonance) 


V  _ { cm  ) 


S.  (cm  mol  cm  ) 


2056.986 

2.56 

X 

10 

6.618 

X 

10 

2057.054 

5.77 

X 

10“^^ 

1.554 

X 

10 

2056.982 

9.35 

X 

10-23 

2.361 

X 

10 

2056.081 

3.13 

X 

10-23 

7.410 

X 

10 

2056.927 

5.81 

X 

10-23 

9.45 

X 

10 

2057.135 

5.89 

X 

10-23 

8.445 

X 

10 

2056.847 

8.94 

X 

10-24 

7.208 

X 

10 

2057.163 

1.65 

X 

10-22 

2.298 

X 

10 

2056.831 

2.08 

X 

10-23 

1.476 

X 

10 

2057.174 

1.09 

X 

10-23 

1.103 

X 

10 

2056.778 

4.44 

X 

10-23 

2.140 

X 

10 

2057.235 

5.37 

X 

10-23 

3.321 

X 

10 

2056.773 

6.74 

X 

10-24 

3.142 

X 

10 

2057.284 

1.49 

X 

10-22 

6.546 

X 

10 

2056.767 

4.23 

X 

10-23 

1.895 

X 

10 

2057.299 

1.40 

X 

10-23 

5.588 

X 

10 

2058.803 

1.40 

X 

10-22 

1.547 

X 

10 

2055.160 

3.78 

X 

10-23 

3.839 

X 

10 

2058.835 

1.74 

X 

10-23 

1.85 

X 

10 

2055.135 

3.43 

X 

10-23 

3.362 

X 

10 

2058.914 

8.20 

X 

10-23 

8.029 

X 

10 

2055.000 

1.03 

X 

10 

8.795 

X 

10 

2058.971 

2.90 

X 

10“^^ 

2.676 

X 

10 

2054.929 

2.76 

X 

10-23 

2.200 

X 

10 

2059.022 

1.54 

X 

10-23 

1.349 

X 

10 

2054.782 

1.27 

X 

10-23 

8.840 

X 

10 

2059.086 

1.70 

X 

10-22 

1.399 

X 

10 

2054.754 

1.01 

X 

10-22 

6.858 

X 

10 

2059.124 

1.00 

X 

10-22 

7.933 

X 

10 

2054.722 

1.56 

X 

10-23 

1.030 

X 

10 

2059.180 

8.11 

X 

10-24 

2.104 

X 

10 

2054.621 

3.84 

X 

10-23 

2.327 

X 

10 

2059.188 

1.97 

X 

10-22 

1.472 

X 

10 

2054.558 

6.76 

X 

10-24 

3.890 

X 

10 

2059.208 

5.90 

X 

10-24 

4.327 

X 

10 

2054.498 

1.26 

X 

10-23 

6.910 

X 

10 

2059.251 

1.51 

X 

10-23 

1.065 

X 

10 

2054.466 

9.00 

X 

10-24 

4.813 

X 

10 

2059.328 

2.07 

X 

10-22 

1.364 

X 

10 

2054.436 

1.77 

X 

10-23 

9.248 

X 

10 

2059.351 

5.92 

X 

10-23 

3.825 

X 

10 

2054.430 

9.97 

X 

10-23 

5.185 

X 

10 

2059.359 

1.36 

X 

10-23 

8.728 

X 

10 

2054.395 

8.40 

X 

10-23 

4.253 

X 

10 

2059.376 

2.14 

X 

10-22 

1.354 

X 

10 

2054.360 

8.53 

X 

10-24 

4.206 

X 

10 

2059.447 

1.87 

X 

10-23 

108 

1.115 

X 

10 
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.J*. 
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2054.295 

6.51 

X 

10-23 

3.059 

X 

10-25 

2059.507 

5.09 

X 

10-2^ 

2.889 

X 

10-25 

2054.292 

1.19 

X 

10-22 

5.580 

X 

10-25 

2089.519 

5.28 

X 

10-23 

2.968 

X 

10-25 

2054.290 

4.12 

X 

10-23 

1.929 

X 

10-25 

2059.523 

8.26 

X 

10-24 

4.268 

X 

10-25 

2054.216 

3.92 

X 

10-23 

1.740 

X 

10-25 

2059.545 

6.04 

X 

10-24 

3.326 

X 

10-25 

2054.171 

2.35 

X 

10-23 

1.010 

X 

10-25 

2059.597 

9.47 

X 

10-23 

5.00 

X 

10-25 

2054.141 

3.93 

X 

10-23 

1.655 

X 

10-25 

2059.691 

3.95 

X 

10-23 

1.943 

X 

10-25 

2054.073 

3.92 

X 

10-23 

1.575 

X 

10-25 

2059.704 

5.73 

X 

10-23 

2.792 

X 

10-25 

2054.030 

7.49 

X 

10-23 

2.924 

X 

10-25 

2059.707 

6.80 

X 

10-23 

3.306 

X 

10-25 

2054.028 

6.64 

X 

10-23 

2.589 

X 

10-25 

2059.715 

3.69 

X 

10-23 

1.783 

X 

10-25 

2053.987 

5.37 

X 

10-24 

2.038 

X 

10-25 

2059.736 

5.23 

X 

10-24 

2.489 

X 

10-25 

2053.944 

1.16 

X 

10-23 

4.480 

X 

10-25 

2053.944 

1.42 

X 

10-23 

5.240 

X 

10-25 

2059.770 

2.22 

X 

10-23 

1.030 

X 

10-25 

2059.879 

1.29 

X 

10-22 

5.539 

X 

10-25 

2053.907 

4.11 

X 

10-23 

1.481 

X 

10-25 

2059.942 

4.10 

X 

10-23 

1.685 

X 

10-25 

5;:'- 


'•'S- 

•I-Xs-; 


2053.859 

2.14 

X 

10 

7.477 

X 

10 

2059.948 

1.37 

X 

10-23 

5.608 

X 

10 

2053.809 

1.87 

X 

10-23 

6.332 

X 

10 

2059.964 

6.03 

X 

10-23 

2.442 

X 

10 

2053.800 

1.47 

X 

10-23 

4.950 

X 

10 

2060.037 

2.97 

X 

10-23 

1.145 

X 

10 

2053.783 

1.07 

X 

10-23 

3.565 

X 

10 

2060.115 

2.13 

X 

10-22 

7.803 

X 

10 

2053.553 

1.01 

X 

10-22 

2.934 

X 

10 

2053.463 

1.08 

X 

10-22 

2.898 

X 

10 

2060.258 

1.64 

X 

10-22 

5.488 

X 

10 

2060.348 

1.78 

X 

10-22 

5.639 

X 

10 

2060.616 

1.98 

X 

10-22 

5.372 

X 

10 

2056.697 

6.86 

X 

10-23 

2.018 

X 

10 

2057.309 

2.33 

X 

10-23 

8.742 

X 

10 

2056.609 

3.98 

X 

10-23 

7.560 

X 

10 

2057.340 

1.64 

X 

10-23 

5.129 

X 

10 

2056.546 

1.21 

X 

10-23 

1.761 

X 

10 

2057.354 

1.19 

X 

10-22 

3.444 

X 

10 

2056.534 

3.13 

X 

10-23 

4.346 

X 

10 

2057,387 

9.34 

X 

10-23 

2.273 

X 

10 

2056.519 

9.67 

X 

10-24 

1.167 

X 

10 

2056.456 

1.47 

X 

10-22 

1.577 

X 

10 

2057.543 

7.54 

X 

10-24 

9.663 

X 

10 

2057.543 

1.20 

X 

10-23 

1.538 

X 

10 

2056.445 

1.51 

X 

10-22 

1.522 

X 

10 

2057.605 

1.57 

X 

10"^" 

1.610 

X 

10 

2056.403 

9.59 

X 

10-2^ 

8.621 

X 

10" 

2057.697 

1.11 

X 

10-23 

8.581 

X 

10" 

2056.352 

1.45 

X 

10-23 

1.115 

X 

lo" 

2057.711 

1.43 

X 

10-22 

1.034 

X 

10" 

2056.217 

1.41 

X 

10-22 

7.536 

X 

10" 

2057.733 

8.94 

X 

10-2^ 

6.180 

X 

10" 

2056.196 

1.58 

X 

10-23 

8.024 

X 

10" 

2057.907 

1.67 

X 

10-23 

7.373 

X 

10" 

2056.125 

1.22 

X 

10-22 

5.260 

X 

10' 

2057.970 

1.59 

X 

10-23 

6.194 

X 

10' 

2055.971 

9.19 

X 

10-2^ 

2.892 

X 

10' 

2057.991 

1.68 

X 

10-“ 

6.265 

X 

10' 

2055.924 

1.44 

X 

10-23 

4.154 

X 

10' 

2058.058 

3.45 

X 

10-23 

1.126 

X 

10’ 

2055.878 

2.03 

X 

10-23 

5.396 

X 

10' 

2058.071 

4.60 

X 

10-23 

1.464 

X 

10 

2055.867 

1.00 

X 

10-22 

2.608 

X 

10 

2058.101 

1.65 

X 

10-22 

4.964 

X 

10 

2055.772 

9.68 

X 

10-2^ 

2.157 

X 

10 

2058.122 

3.53 

X 

10-23 

1.022 

X 

10 

2055.751 

1.42 

X 

10-“ 

3.061 

X 

10 

2058.152 

7.48 

X 

10-2^ 

2.052 

X 

10 

2055.733 

4.60 

X 

10-23 

9.642 

X 

10 

2058.154 

8.08 

X 

10-23 

2.209 

X 

10 

2035.718 

1.31 

X 

10-22 

2.683 

X 

10 

2058.243 

4.93 

X  10-23 

1.1587X  10-2^ 

2055.683 

8.50 

X  10-23 

1.651 

X  10-24 

2058.312 

2.16 

X  10-23 

4.548 

X  10-23 

2055.632 

2.72 

X  10-23 

4.905 

X  10-23 

2058.335 

6.88 

in-24 

X  10 

1.398 

X  10-23 

2055.572 

5.88 

in~24 

X  10 

9.748 

10-26 
X  10 

2058.369 

9.30 

in“23 

X  10 

1.796 

10-24 
X  10 

2055.489 

4.88 

ift-23 

X  10 

7.240 

X  10-23 

2058.380 

5.19 

ift-23 

X  10 

9.862 

X  10-23 

2055.484 

7.70 

irt-23 

X  10 

1.135 

X  10-24 

2058.390 

2.19 

in-23 

X  10 

4.101 

X  10-23 

2055.432 

2.35 

X  10-23 

3.242 

X  10-23 

2058.517 

3.59 

X  10-23 

5.629 

X  10-23 

2055.384 

6.03 

X  10-2^ 

7.839 

10-26 
X  10 

2058.521 

2.77 

X  10-23 

9.320 

10-24 
X  10 

2055.355 

3.87 

X  10-23 

4.858 

X  10-23 

2058.566 

1.23 

X  10-2^ 

1.808 

X  10-23 

2055.333 

1.98 

X  10-23 

2.42 

X  10-23 

2058.576 

5.24 

X  10-23 

7.61 

X  10-23 

2055.321 

1.15 

X  10-23 

1.386 

10-25 
X  10 

2058.596 

1.94 

X  10-22 

2.744 

10-24 

X  10 

2055.319 

6.36 

X  10-23 

7.650 

10-25 
X  10 

2058.603 

6.78 

X  10-23 

9.507 

X  10-23 

2055.276 

3.60 

X  10-23 

4.120 

X  10-23 

2058.627 

1.86 

X  10-22 

2.531 

X  10-24 

2055.246 

1.02 

X  10-22 

1.128 

X  10-24 
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I 


m 


K'-i" 


2058.724 

6.24 

X 

10 

7.550 

X 

10 

2055.225 

6.80 

X 

10-23 

7.348 

X 

10 

2058.728 

2.65 

X 

10-23 

3.191 

X 

10 

2055.220 

1.78 

X 

10-23 

1.913 

X 

10 

2058.744 

1.86 

X 

10-23 

2.199 

X 

10 

2055.162 

1.17 

X 

10-22 

1.180 

X 

10 

2060.677 

1.07 

X 

10-22 

2.807 

X 

10 

2060.798 

2.41 

X 

10-22 

5.923 

X 

10 

2060.869 

1.05 

X 

10-22 

2.486 

X 

10 

2060.991 

2.30 

X 

10-22 

5.116 

X 

10 

2061.026 

1.11 

X 

-22 

10 

2.426 

X 

10 

2061.324 

1.10 

X 

10-22 

2.083 

X 

10 

2061.355 

1.41 

X 

10-22 

2.631 

X 

10 

2061.488 

20.58X 

10-22 

4.532 

X 

10 

2061.683 

1.92 

X 

10-22 

3.097 

X 

10 

2061.813 

1.88 

X 

10-22 

2.870 

X 

10 

2061.989 

2.33 

X 

10-22 

3.301 

X 

10 

2062.109 

2.71 

X 

10-22 

3.669 

X 

10 

2062.158 

2.87 

X 

10-22 

3.812 

X 

10 

Total  <y 

-  1.37 

X 

10 

“fc?  rv« 
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Appendix  H 


CO2  Extinction  Coefficient  ( "Off "-Resonance) 


a —  (cm  ) 


(cm"^) 


X  10  ^(]an"^) 


1.361 

X 

10 

8.405 

X 

10-^ 

1.144 

7.200 

X 

10-22 

7.626 

X 

10^3 

.605 

4.94 

X 

10-22 
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H2O  Extinction 


Coefficient  ( "Off "-Resonance) 
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Appendix  L 

Hot  C0»  Absorption  Cross  Section  ( "Off-Resonance) 
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